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CHAPTER  1 
INTRODUCTION  AND   LITERATURE  REVIEW 

1.1.  Motivation  for  the  study.  An  important  phase  of  dairy 
herd  management  has  to  do  with  elimination  of  low  producing 
cows.  Because  of  impracticality,  herds  cannot  begin  with  all 
cows  being  at  exactly  the  same  stage  of  lactation  at  the  same 
time  nor  will  they  all  be  the  same  age.  Many  times  it  is  also 
economically  impractical  to  wait  until  all  cows  have  completed 

a  lactation  before  making  culling  decisions.  For  these  reasons, 
it  is  necessary  to  find  some  way  of  converting  the  performance 
of  all  cows  in  a  herd  to  a  standard  basis  with  respect  to  milk 
and  butterfat  production.  An  obvious  procedure  is  to  convert 
the  production  of  each  cow  to  expected  production  given  the 
same  environment,  level  of  maturity,  and  a  completed  lactation 
of  305  days.   One  could  then  compare  all  cows  on  the  same  basis. 

1.2.  Literature  concerning  factors  for  converting  records  to 
505  days. 

1.2.1.  Preliminaries.  Since  converting  cows  to  a  standard 
basis  for  comparison  is  not  a  new  problem,  and  since  there  are 
other  applications  of  the  extension  of  partial  production 
records  to  305  days  such  as  sire  proving,  evaluating  sold  and 
dead  cows,  etc.,  tables  of  factors  for  converting  milk  and 
fat  production  to  305  days  from  a  partial  record  have  already 
been  prepared,  and  the  effects  of  certain  environmental 
variables  upon  the  shape  of  the  lactation  curve  as  expressed 
by  these  factors  have  been  studied. 


1.2.2.  Types  of  factors.  Environmental  variables  have  been 
studied  in  terms  of  three  types  of  factors.  The  first  type 
consists  of  those  for  finding  a  quick  estimate  of  the  producing 
ability  of  a  cow  when  full  information  is  lacking.  Such 
factors  were  determined  from  which  a  one-day  test  could  be  used 
to  estimate  producing  ability.  Adjustments  were  presented  by 
Cannon  et  al.  [1942],  Madden  et  al.  [1956,  1959],  Lamb  and 
McGilliard  [1960a],  and  Van  Vleck  and  Henderson  [1961b].  The 
second  type  is  used  for  estimating  total  production  when 
several  months  records  are  available.  Van  Vleck  and  Henderson 
[1961b,  e]  presented  a  set  of  regression  coefficients  to  be 
used  when  m  sequential  monthly  records  were  available,  either 
at  the  beginning  or  the  end  of  the  lactation.  The  third  type 
consists  of  those  based  on  cumulative  production.   Cumulative 
factors  were  presented  by  Cannon  et  al.  [1942],  Madden  et  al. 
[1956,  1959],  lamb  and  McGilliard  [1960a],  and  Van  Vleck  and 
Henderson  [1961b,  c,  e].  This  third  type  is.  of  most  interest 
to  the  present  study. 

1-2. 3-  Variables  studied.  The  cumulative  factors  mentioned 
above  were  studied  in  terms  of  breed,  age,  lactation  number, 
milking  frequency,  level  of  production,  herd,  and  season,  though 
not  all  variables  were  examined  in  all  the  above  studies.   Of 
all  variables  considered,  age  was  most  important.  Factors 
could  be  classified  into  two  groups,  i.e.,  those  for  cows  less 
than  3  years  old,  and  those  for  cows  greater  than  or  equal  to 
3  years  of  age.   It  was  found  that  even  when  records  were 
adjusted  for  age  before  determining  the  coefficients,  all  age 


effects  were  still  not  removed,  and  thus,  two  corrections  were 
needed,  one  to  adjust  for  the  effect  of  age  on  the  shape  of 
the  lactation  curve,  and  the  other  to  adjust  for  differences  in 
level  of  production  between  older  and  younger  cows.   [See 
Madden  et  al. ,  1955,  1959,  and  Lamb  and  McGilliard,  1960a.] 
Van  Vleck  and  Henderson  [1961b,  c]  also  found  it  necessary  to 
adjust  for  age.  Little  difference  was  found  by  Lamb  and 
McGilliard  [1960a]  between  the  effect  of  age  and  that  of 
lactation  number. 

Breed  was  found  to  be  significant  by  Lamb  and  McGilliard 
[1960a],  while  Cannon  et  al.  [194-2]  found  that  because  of  the 
similarity  among  regression  coefficients  for  the  various  breeds 
studied,  the  coefficients  could  be  combined.  Madden  et  al. 
[1955,  1956,  1959]  and  Van  Vleck  and  Henderson  [1961b,  c,  e] 
studied  only  Holsteins. 

Lamb  and  McGilliard  [1960a]  and  Van  Vleck  and  Henderson 
[1961b,  c]  found  season  to  be  significant,  and  Lamb  and 
McGilliard  found  a  breed  by  season  interaction,  with  Holsteins 
being  least  influenced  by  season  of  freshening. 

Madden  et  al.  [1959]  found  milking  frequency  and  level 
of  production  to  be  nonsignificant.  Lamb  and  McGilliard 
[1960a]  found  herd  differences  to  be  non-significant,  while 
Van  Vleck  and  Henderson  [1961b,  e]  presented  their  factors 
both  on  a  wi thin-herd  basis  and  ignoring  herd  effects.   In  a 
later  study,  Van  Vleck  and  Henderson  [1961e]  found  that  in 
estimating  the  regression  ignoring  herds  [on  a  total  regression 
rather  than  on  an  intra-herd  regression  basis]  the  inaccuracy 


due  to  ignoring  herd  differences  was  greatly  outweighed  by  the 
difficulty  in  obtaining  the  herd  means  and  the  inaccuracy  of 
their  estimates. 

Most  tables  of  estimates  for  extending  records  are  based 
on  either  milk  or  fat  alone  but  are  used  to  predict  both.   It 
was  found  by  Lamb  and  McGilliard  [1960a]  that  there  are 
definite  differences  between  milk  and  fat  adjustments;  with 
the  exception  of  the  first  month  for  Holsteins  and  Brown  Swiss, 
cumulative  factors  for  milk  will  underestimate  butterfat 
production  for  the  lactation.  They  also  found  that  season  of 
freshening  was  important  for  milk,  but  not  for  fat.  Van  Vleck 
and  Henderson  [1961b]  found  coefficients  for  milk  and  fat  to 
be  distinctly  different.  Madden  et  al.  [1955]  presented 
slightly  different  factors  for  milk  and  fat.  Van  Vleck  and 
Henderson  [1961c]  raised  the  question  whether  milk  and  fat 
should  both  be  extended,  as  this  could  change  the  per  cent  fat 
for  the  total  lactation.  The  systems  now  used  do  not  affect 
the  fat  percentage  since  they  are  based  solely  on  milk  or  on 
fat. 

1.2.4.  Procedures  Used.   In  the  papers  mentioned  above  which 
attempt  to  evaluate  effects  of  certain  variables  on  milk  and 
fat  production,  two  approaches  are  used:   [1]  a  ratio  factor 
is  constructed  based  upon  either  the  total  milk  produced  to 
date  divided  into  the  total  milk  produced  in  305  days,  the 
total  milk  produced  on  x  consecutive  test  days  divided  into 
the  total  milk  produced  on  the  10  consecutive  test  days,  or  the 
sum  of  the  reciprocals  of  the  non-cumulative  factors,  the  sum 


then  being  reciprocated,  for  each  x  periods.  [By  non- cumulative 
is  meant  the  per  cent  of  the  total  milk  produced  each  month  of 
the  lactation,  or  on  each  test  day  of  the  lactation.]  This 
will  give  a  cumulative  ratio  factor.   [2]  a  regression 
coefficient,  expressed  either  in  terms  of  a  regression  equation 
or  a  single  coefficient,  is  found  from  the  regression  of  total 
milk  produced  in  305  days  on  total  milk  produced  to  date.   In 
each  of  these  methods  the  factor  is  adjusted  in  some  way  for 
such  variables  as  age,  season  of  freshening,  parity,  and  breed. 

The  comparative  accuracy  of  the  above  two  approaches  has 
been  discussed  by  Madden  et.  al.  [1959].  They  found  that  the 
ratio  method  may  underestimate  the  total  production  of  low 
producing  cows  and  over— estimate  total  production  of  high 
producing  cows  since  the  ratio  method  adjusts  only  for  the 
incompleteness  of  the  lactation  and  does  not  take  into  account 
the  incomplete  repeatability  of  the  parts  of  the  lactation, 
which  is  a  part  of  the  regression  method.  This  difference  is 
largest  during  the  early  months  of  a  lactation.  The  regression 
method,  in  which  the  cumulative  part  of  the  lactation  is  multi- 
plied by  the  regression  [b]  of  whole  on  part  and  added  to  an 
appropriate  constant  [a]  to  estimate  the  total  production,  not 
only  adjusts  for  incompleteness,  but  also  for  unidentified 
sources  of  variation  which  make  the  part  greater  than  or  less 
than  the  average.  The  total  estimated  by  regression  varies  less 
than  the  actual  total.  Although  the  variance  of  total  production 
estimated  by  the  ratio  method  is  greater  than  that  from  the 
estimate  found  by  regression,  it  has  variation  similar  to 


actual  total  production.  For  this  reason,  if  these  extended 
records  are  to  be  used  for  culling  purposes  or  sire  proving, 
ratio  techniques  are  to  be  preferred.   Since  records  extended 
by  regression  differ  less  than  actual,  this  tendency  to  group 
the  records  around  the  mean  may  make  selection  decisions  more 
difficult. 

1.2.5.  Correlations  of  individual  month  records  with  complete 
record.  The  part  and  the  whole  are  always  correlated;  Van  Vleck 
and  Henderson  [1961b]  give  the  correlation  of  the  individual 
month's  record  with  the  complete  record  as  0.57,  O.75,  0.81, 
0.85,  0.85,  0.85,  0.85,  O.78,  0.66,  and  O.53  for  months 
1,...,  10,  respectively,  on  a  within-herd  basis.  The  correlations 
for  the  cumulative  months  with  the  complete  record  [on  a  within- 
herd  basis]  are  0.57,  O.75,  0.82,  0.87,  0.90,  0.93,  0.95,  0.97, 
and  0.99  for  months  1,...,  9,  respectively.  When  herd  is  not 
considered  [see  Van  Vleck  and  Henderson  (1961e)]  the  following 
correlations  are  obtained  for  the  individual  months  0.67,  0.82, 
0.86,  0.89,  0.90,  0.89,  0.87,  0.81,  0.68,  and  0.52  for  months 
1,...,  10,  respectively.  The  correlations  for  the  cumulative 
months  records  are  0.67,  0.82,  0.87,  0.91,  0.93,  0.95,  0.97,  and 
0.99  for  months  1,...,  9,  respectively. 

1.3.   Literature  concerning  factors  for  converting;  records  to 
mature  equivalent. 

I.3.I.  Preliminaries.   If  one  desires  to  convert  all  cows  in  a 
herd  to  the  same  basis,  when  the  cows  are  at  different  stages 
of  lactation,  by  extending  their  production  to  305  days  and 
converting  them  to  mature  equivalent,  one  must  not  only  be 


concerned  with  extrapolation,  but  also  with  age  correction 
adjustments. 

Milk  production  increases  with  age  at  an  ever-decreasing 
rate  until  maximum  production  is  reached  at  6  to  8  years;  it 
then  declines  with  advancing  age.  The  regression  of  production 
on  age  is  distinctly  curvilinear.  Lush  and  Shrode  [1950]  have 
given  the  theory  and  problems  of  age  adjustments  and  their 
calculation.  They  described  two  methods  and  the  bias  associated 
with  each.  When  the  age  corrections  are  calculated  by  Method  A, 
averaging  all  records  made  at  each  age,  concurrent  selection 
introduces  a  bias.   If  the  factors  are  calculated  by  Method  B, 
comparing  only  records  made  by  the  same  cows  at  two  successive 
ages,  concurrent  selection  introduces  a  bias  in  the  opposite 
direction.  The  true  age-change  would  then  lie  somewhere  between 
the  apparent  changes  given  by  the  two  methods. 

Lush  and  Shrode  also  bring  out  the  fact  that  a  bias  is 
introduced  if  the  average  inherent  productivity  of  the  dairy 
population  is  increasing.  For  instance,  at  any  given  date  the 
averages  of  the  older  cows  do  not  yet  include  the  records  from 
cows  born  in  the  most  recent  years  when  the  production  was 
actually  higher.  Henderson  et  al.  [19591  presented  techniques 
to  estimate  this  genetic  time  trend  with  age. 

1.3.2.  Studies  of  the  effects  of  age  on  production.  The  effects 
of  age  on  production  have  been  studied  by  workers  other  than 
Lush  and  Shrode  [1950].  Searle  and  Henderson  [1959,  I960], 
Searle  [I960,  1961a]  and  Van  Vleck  and  Henderson  [1961c]  have 


presented  various  approaches  and  techniques  to  the  calculation 
of  adjustments  for  the  effects  of  age  on  305-day  production. 
I.3.3.  Which  factors  should  he  used.   Of  all  the  adjustments 
in  present  use,  the  standard  DHIA  factors  given  by  Kendrick 
[19552  are  probahly  the  most  widely  used.  Miller's  study 
[Killer,  1964]  would  indicate  that  in  the  absence  of  further 
research,  these  factors  would  be  the  wisest  to  use.  Thus  no 
attempt  was  made  to  estimate  age-correction  factors  in  the 
present  study.  Other  reasons  will  be  given  in  Section  1.5- 
1.4.  Literature  concerning  statistical  considerations. 

1.4.1.  Preliminaries.   In  the  statistical  analyses  of  variables 
which  affect  the  shape  of  the  lactation  curve ,  there  are  several 
conditions  which  cause  problems  with  respect  to  the  use  of 
standard  analysis  of  variance  and  regression  techniques. 

1.4.2.  Non-orthogonal  data.  When  dealing  with  analyses  in 
which  the  numbers  in  each  subclass  are  unequal,  the  problem 
becomes  difficult,  because  standard  analyses  of  variance 
techniques  assume  that  degrees  of  freedom  on  which  the  variances 
for  each  of  the  effects  being  analyzed  are  based  are  orthogonal 
to  [or  independent  of]  each  other.   When  unequal  subclass 
frequencies  are  encountered,  orthogonality  cannot  be  guaranteed. 
Several  papers  have  been  written  proposing  either  exact  or 
approximate  solutions  to  this  problem.   [See  Yates,  1934; 
Snedecor,  1934;  Harvey,  I960;  Gosslee  and  Lucas,  1965;  Wakefield, 
1965;  and  Mielke  and  McHugh,  1965]-  Several  of  these  techniques 
are  listed  and  discussed  in  various  statistical  books,  for 
example,  Snedecor  [1956],  Goulden  [1952],  Brownlee  [I960], 


Graybill  [1961],  and  Dixon  [19651- 

Among  techniques  presented  by  Yates  [19540  was  a  simple 
approximate  method  which  consisted  of  ignoring  differences 
between  subclass  numbers  and  performing  an  analysis  of  subclass 
means,  assuming  the  variances  of  the  means  to  be  equal  and 
considering  the  means  as  one  observation  per  cell.   The 
analysis  is  then  performed  using  techniques  for  equal  sub- 
classes. This  approximation  is  only  useful,  he  pointed  out, 
if  the  subclass  numbers  do  not  differ  greatly.  Recently, 
Gosslee  and  Lucas  [1965]  concluded,  concerning  the  question  of 
level  of  significance  in  the  method  of  unweighted  squares  of 
means,  that  the  effects  of  unequal  variance  among  the  cell 
means  have  only  a  moderate  effect. 

Estimates  of  components  of  variance  are  also  affected  by 
unequal  subclass  numbers.  Methods  of  dealing  with  this  problem 
have  been  presented  by  Henderson  [1953],  Harvey  [I960],  and 
Bush  and  Anderson  [1963].  However  in  the  case  of  unweighted 
squares  of  means  [Yates,  1934,  and  Gosslee  and  Lucas,  1965], 
it  is  necessary  only  to  solve  the  expected  mean  squares  in  the 
usual  manner  for  orthogonal  data. 

1.4.3-  Multiple  comparisons.  Use  was  made  of  Duncan's 
Multiple  Range  Test  to  test  for  differences  among  means  in 
several  of  the  papers  discussed  in  Sections  1.2  and  1.4.2. 
Although  Duncan's  test  decreases  the  probability  of  a  Type  II 
error,  the  probability  of  a  Type  I  error  inflates  until  no 
estimate  of  it  even  exists.  For  this  reason,  tests  for 
which  the  probability  of  the  Type  I  error  remains  constant 
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and  known,  such  as  Tukey's  hsd  or  Sheffe's  test,  are  to  be  pre- 
ferred, [see  Federer,  1955,  Brownlee,  I960,  and  Sheffe',  1959]. 
1.4-.4.  The  comparison  of  growth  curves.   Comparing  several 
lactation  curves  is  similar  to  the  problem  of  comparing  growth 
curves.  Rao  [1958]  defines  this  problem  in  terms  of  comparing 
the  characteristics  of  growth  under  different  conditions  such 
as  diet,  environment,  etc.  The  difficulty  in  comparing  growth 
curves  based  on  values  at  a  number  of  points  along  the  curve  can 
be  overcome  by  reducing  the  data  to  the  lowest  possible  dimensions 
without  sacrificing  the  essential  information.  This  can  lead 
to  a  more  efficient  procedure.  As  an  example  of  this,  he  cites 
the  fitting  of  a  second  degree  polynomial  to  each  growth  curve 
classified  by  the  factors  to  be  studied.  The  coefficients  of 
the  linear  and  quadratic  terms  can  then  be  taken  to  represent 
the  salient  features  of  the  growth,  and  thus  a  large  number  of 
observations  can  be  replaced  by  two  coefficients.  A  small  number 
of  parameters  to  use  in  the  determination  of  differences  between 
groups  of  growth  curves  can  thus  be  found. 

1.4. 5.  The  distribution  of  the  ratio  of  two  normally  distributed 
random  variables.  The  lactation  estimation  curve  is  generally 
expressed  as  a  series  of  10  ratios  between  the  total  milk 
produced  to  date  and  the  total  305-day  milk  production,  where 
these  two  values  become  more  and  more  highly  correlated  as  the 
ratios  are  taken  from  the  first  to  the  tenth  period,  being 
exactly  1.00  for  the  tenth  period.  These  ratios  are  then  ratios 
of  two  random  variables  which  are  assumed  to  be  normal.  The 
distribution  of  the  ratio  of  two  normally  distributed  random 
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variables  is  not  normal.   Geary  [1930],  Fieller  [1932], 

Merrill  [1928],  and  Marsaglia  [1965]  have  considered  this 

problem  in  general,  and  in  its  specific  applications,  such,  as 

bioassay  and  sampling  theory,  it  has  been  discussed  by  Fieller 

[1940],  Finney  [1964],  Hansen,  et  al.  [1956],  and  Deming 

[1950,  I960]. 

Although  the  distribution  of  z  =  y/x  is  not  normal, 

there  exists  a  function  of  z  which  is  approximately  normal. 

If  x,  y  are  distributed  normally  with  correlation  r  and 

with  means  0  and  standard  deviations  a  and  0,  respectively, 

then  if   z  ■     «g  ,   where   a  and  b  are  non- negative 

constants,  then  t  =  az  is  distributed 

v  a2  z2  +  2ra|3z  +  g2 
H[0,1]   if  a  +  x  is  unlikely  to  assume  negative  values.  This 

condition  is  amply  satisfied  if  a  >  30c,  i.e.,  if  the  coefficient 

of  variation  of  a  +  x  is  not  greater  than  1/3-  Thus,  normal 

theory  can  be  used  when  ratios  make  up  the  observations  and  the 

x  and  y  values  are  positive  and  not  close  to  zero  and  are 

large  with  respect  to  their  standard  deviations,  provided  the 

above  transformation  is  made. 

1.5>  The  purpose  of  this  thesis.  The  purpose  of  this  thesis  is 

to  examine  the  effects  of  certain  variables  upon  the  shape  of 

the  lactation  estimation  curve,  using  a  different  approach  to 

the  problem  than  those  which  have  been  heretofore  presented, 

taking  into  account  the  findings  presented  in  the  various 

sections  of  this  introduction.  The  data,  from  herds  in  the 

Western  states — data  which  have  not  previously  been  available, 

will  be  used  to  provide  factors  based  upon  the  significant 
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variables  for  use  in  estimating  total  from  part  lactation 
production  for  both  milk  and  fat. 

Because  of  the  nature  of  the  data  available,  it  was  not 
practical  to  construct  factors  for  converting  cow  records  to 
mature  equivalent.  For  this  reason,  and  because  of  the 
availability  of  the  age  correction  factors  already  in  existence, 
it  was  decided  not  to  consider  them  in  the  present  study. 


CHAPTER  2 
METHODS  AND  MATERIALS 

2.1.  Materials.  Data  used  in  this  study  came  from  the  files 
of  the  DHI  Computing  Service,  Provo,  Utah,  and  comprised  all 
records  tabulated  by  the  computing  service  from  November,  1962, 
through  March,  1965,  excluding  records  for  the  state  of  California. 
1,400,800  monthly  cow  production  totals  were  recorded  on  magnetic 
tape  and  sorted  into  sequence  by  month  of  production  within  cow. 
The  records  were  then  edited  to  satisfy  the  following  conditions. 

1.  A  lactation  consisting  of  305  to  350  days.   [Note  that 
when  a  lactation  exceeded  305  days,  the  production  was  converted 
back  to  305  days.]  There  is  some  question  whether  to  use  short 
lactation  records,  i.e.,  completed  lactations  of  cows  which  have 
been  in  milk  more  than, say,  270  days  but  less  than  305,  because 
a  large  number  of  cows  complete  their  lactations  before  reaching 
305  days. 

2.  No  missing  test  days  during  the  lactation.  Sometimes, 
if  a  cow  is  purchased  fresh,  the  first  part  of  her  lactation 
will  be  missing.  Although  it  was  highly  unlikely  because  of 
the  way  in  which  the  records  were  assembled,  test  dates  were 
tested  for  consecutiveness  to  eliminate  lactations  with  non- 
consecutive  dates. 

3.  The  first  month's  test  being  used  to  estimate  less 
than  or  equal  to  50  days'  production.  This  restriction  was 
used  to  avoid  overlapping  of  effects  from  the  first  and  second 
period's  production,  as  this  could  mask  or  confound  the  effects 
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of  number  of  days  into  the  lactation  on  the  lactation  estimation 
curve.  At  the  DHI  Computing  Service,  as  soon  as  a  cow  exceeds 
320  days  in  milk  without  completing  her  lactation,  i.e. 
approximately  one  month  after  reaching  305  days,  the  lactation 
is  automatically  flagged  and  a  305  day  record  computed.  Because 
of  the  choice  of  50  days  as  the  cut-off  for  the  first  period, 
350  days  becomes  the  maximum  possible  number  of  days  in  the 
flagged  lactation.  As  was  mentioned  before,  such  records  were 
then  converted  back  to  305  days. 

4--  The  latest  record  being  used  if  duplicate  records  appeared 
for  any  given  month.  It  was  felt  that  if  such  duplicates  appeared, 
the  latest  record  would  be  the  most  likely  to  be  correct. 

5-  Two-time  milking.  There  were  not  enough  three  times 
per  day  milking  records  to  make  their  consideration  worthwhile. 

6.   A  lactation  number  reported.   In  these  records,  the 
age  of  the  cow  was  not  available.  Since  age  is  one  of  the  vari- 
ables to  be  studied  [see  sections  1.2.2,  and  2.2.1],  a  record  was 
worthless  unless  some  estimate  of  age  was  available,  thus  lacta- 
tion number  was  necessary. 

7-  Wo  correction  appearing  anywhere  in  the  lactation. 
Often  in  the  reporting  of  information,  mistakes  and  inconsisten- 
cies occur.  When  these  mistakes  are  found,  a  correction  card  is 
placed  in  the  file  to  adjust  the  record  for  the  cumulative  pro- 
duction.  A  correction  card  is  also  used  to  estimate  missing 
production,  correct  for  sickness,  etc.  Because  these  corrections 
indicate  incorrect  data  in  the  early  part  of  the  lactation  which 
could  affect  the  shape  of  the  lactation  estimation  curve,  records 
containing  them  were  eliminated. 
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8.  Consistency  of  data.  Whenever  data  were  obviously 
incorrect  or  inconsistent,  i.e.,  milk  weights  that  were  too 
large  [99999  in  that  portion  of  the  record,  etc.],  or  when 
data  needed  for  the  analysis  were  missing,  such  lactations 
were  eliminated. 

It  should  he  noted  with  respect  to  these  data  that  not 
only  was  the  age  of  the  cow  not  available,  hut  also,  it  was 
not  possible  to  determine  breed  from  the  data. 

For  the  analysis  [which  will  be  referred  to  later]  of 
Days  Dry  and  Days  Open,  a  further  edit  was  performed.   In 
addition  to  the  above  conditions,  those  data  used  in  this 
analysis  were  also  required  to  have  a  non-zero  number  of  days 
carried  calf  and  be  second  lactations  or  greater. 

There  were  18,541  records  satisfying  the  above  8  conditions, 
and  13,023  satisfying  also  the  conditions  mentioned  in  the 
preceding  paragraph. 

Examples  of  the  form  of  the  data  before  and  after  editing 
appear  in  Appendix  I. 
2.2.   Methods. 

2.2.1.  Variables  to  be  studied.   It  was  decided  to  compare 
lactation  estimation  curves  based  on  the  edited  data  with 
respect  to  the  following  variables. 

1.  Geographic  Area.   It  was  decided  to  examine  the  effects 
of  different  groupings  of  the  states  represented  in  these  data. 

2.  Average  per  cent  fat.  Because  breed  was  not  available 
to  be  used  as  a  variable  in  this  study,  average  per  cent  fat 
was  substituted.  Average  per  cent  fat  is  defined  as  total 
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pounds  of  milk  divided  into  total  pounds  of  fat,  giving  the 
weighted  average  over  the  entire  lactation.  An  interesting 
consideration  of  the  use  of  this  variable  rather  than  breed 
comes  from  the  question  of  what  the  real  breed  differences  are 
in  terms  of  milk  and  fat  production.  For  example,  what  is  the 
difference  between  a  Guernsey  that  tests  3-65%  and  a  Holstein 
that  tests  3.65%  fat?  Or  a  Holstein  that  tests  4.75%  and  a 
Guernsey  that  tests  4.75%  fat? 

3-  Lactation  number.   In  all  studies  referred  to  in 
Section  1.2  in  which  they  were  considered,  both  age  and  lacta- 
tion number  were  found  to  be  significant.   Lamb  and  McGilliard 
[1959]  found  that  except  in  cases  where  age  and  lactation  number 
did  not  coincide,  lactation  number  was  actually  a  better  means 
of  estimation. 

4.  Season  of  Freshening.   In  several  of  the  studies 
discussed  in  Section  1.2,  season  was  found  to  have  an  effect  on 
the  shape  of  the  lactation  estimation  curve. 

5-  Previous  days  dry.  As  Smith  [1959-2  and  Johansson 
[19613  point  out,  previous  days  dry  ont  only  have  an  effect  on 
the  current  lactation,  but  also  on  the  succeeding  one. 

6.  Days  open.  Days  open  is  defined  to  be  the  difference 
between  days  carried  calf  and  days  in  milk  plus  three  days.   In 
a  discussion  of  physiology  of  lactation,  both  Smith  [1959]  and 
Johansson  [1961]  mentioned  that  this  variable  was  found  to  have 
an  effect  on  the  lactation. 

7-   Level  of  milk  production.   In  all  factors  in  use  at 
the  present  time,  level  of  milk  production  is  not  taken  into 
account.  Obviously,  it  is  impossible  to  consider  factors 
which  take  this  into  account  in  terms  of  total  milk  production 
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for  a  given  lactation,  unless  one  uses  level  of  production  for 
the  preceding  lactation,  level  of  production  for  the  first  month, 
or  the  herd  average  as  an  indication  of  level  of  production  of 
the  current  one. 

8.   Level  of  fat  production.   It  was  decided  to  examine 
also  the  level  of  fat  production  to  see  if  the  make-up  of  the 
production  cycle  was  affected  by  differences  in  level  of 
production.   It  is  known  that  the  per  cent  fat  test  decreases 
as  the  milk  production  increases,  but  not  what  happens  to  the 
shape  of  the  curve  as  milk  or  fat  production  increases. 
2.2.2.   Levels  and  coding  of  the  variables  to  be  used.   To 
determine  the  levels  of  the  above  factors  to  be  used,  several 
things  were  taken  into  account,  the  most  important  being  the 
frequencies  of  numbers  of  cows  in  each  grouping  and  the 
physiological  and  environmental  effects  involved. 

The  number  of  records  from  each  state  are  given  in  the 
following  table. 

Table  1 
The  number  of  records  from  each  state  used  in  the  study. 
State  Number  of  Records 


North  Dakota  (42) 

19 

Iowa  (45) 

28 

Nebraska  (47) 
Montana  (81) 

467 

718 

Idaho  (82) 

553 

Wyoming  (83) 

775 

Colorado  (84) 

4 

.715 

New  Mexico  (85) 

1 

,154 

Arizona  (86) 

1 

,440 

Utah  (87) 

a 

,210 

Nevada  (88) 

503 

18,582 
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When  the  analyses  were  begun,  certain  data  were  found  to 
"be  incorrect  and  had  not  been  detected  in  the  original  edit. 
When  these  were  eliminated,  the  number  of  records  was  reduced 
from  18,582  to  18,541.  The  states  were  combined  such  that 
North  Dakota,  Iowa,  Nebraska,  Montana,  Idaho,  and  Wyoming 
were  considered  Area  2;  Colorado  was  considered  Area  3;  New 
Mexico  and  Arizona  were  considered  Area  4;  and  Utah  and  Nevada 
were  considered  Area  5>  These  groupings  were  made  on  the  basis 
of  the  above  frequency  table  [Table  1],  climatic  similarities, 
and  the  necessary  combinations  to  eliminate  missing  cells  when 
the  other  factors  in  the  analysis  were  taken  into  account. 

The  levels  of  average  per  cent  fat  were  determined  from  an 

examination  of  the  average  per  cent  fat  for  the  different 

breeds,  and  also  from  the  mean  and  variance  of  the  frequency 

distribution  for  Utah.   The  mean  was  3.8%  and  the  standard 

deviation  was  0.6%.   When  this  was  done,  the  levels  were 

considered  to  be  those  shown  in  Table  2. 

Table  2 

The  data  classes  of  per  cent  fat  and  the  codes  thereof 
[preliminary] . 

Per  Cent  Fat  Code 

2.8%  and  under  1 

2.9%  to  3-4%  2 

3.5%  to  4.0%  3 

4.1%  to  4.6%  4 

4.7%  to  5-2%  5 

5.3%  to  5.8%  6 

5.9%  and  over  7 

When  these  levels  were  examined  for  cell  frequencies,  it  was 
found  necessary  to  make  the  following  combinations  and  to 
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analyze  the  data  from  Geographic  Areas  3  and  5  Cas  presented 
in  Tables  3  &nd.   4]. 

Table  3 

The  data  classes  of  per  cent  fat  and  the  codes  thereof  for 
Geographic  Areas  2  and  4. 

Per  Cent  Fat  Code 

3.4%  and  under  2 

3.5%  to  4.0%  3 

4.1%  and  over  4 


Table  4 

The  data  classes  of  per  cent  fat  and  the  codes  thereof  for 
Geographic  Areas  3  a*1"!  5> 

Per  Cent  Fat  Code 

3.4-%  and  under  2 

3.5%  to  4.0%  3 

4.1%  to  4.6%  4 

4.7%  to  5.2%  5 

5. 3%  and  over  6 


These  levels  correspond  roughly  to  low  fat  ^producing  Holsteins 
[23,  mostly  Holsteins  and  Milking  Shorthorns  with  some  Ayrshires 
and  Brown  Swiss  C 33 ,  and  a  few  Holsteins  and  Milking  Shorthorns, 
the  remainder  being  Ayrshires,  Brown  Swiss,  Guernseys,  and 
Jerseys  [4]  for  Geographic  Areas  2  and  4,  and  to  low  fat 
producing  Holsteins  [2],  mostly  Holsteins  and  Milking  Shorthorns, 
with  some  Ayrshires  and  Brown  Swiss  [3],  half  of  the  Ayrshires 
and  Brown  Swiss,  and  half  of  the  Guernseys  plus  a  few  Holsteins 
and  Milking  Shorthorns  [4],  Guernseys  and  half  Jerseys  [5],  and 
predominantly  Jerseys  with  a  few  Guernseys  [6]  for  Geographic 
Areas  3  and  5- 

Although  there  has  been  strong  evidence  that,  in  general, 
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age  can  be  divided  on  "less  than  3  years"  and  "3  years  and 

over"  [see  Section  1.2],  it  was  decided  to  look  at  age  [based 

on  lactation  number]  in  terms  of  the  six  levels  given  in  Table 

5. 

Table  5 

The  data  classes  of  lactation  number  and  the  codes  thereof. 

Lactation  Number               Code 

First                        1 
Second                         2 
Third                          3 
Fourth                         4 
Fifth                          5 
Sixth  and  over                 6 

It  was  decided  to  group  season  of  freshening  in  the  manner 

indicated  in  the  following  table. 

Table  6 

The  data  classes  of  season  and  the  codes  thereof. 

Season                        Gtifln 

December,  January 

and  February                 1 

March,  April,  and 

May                         2 

June,  July,  and 

August                      3 

September,  October 

and  November                4 

Previous  days  dry  were  divided  into  six  groups.   Johansson 

[1961]  stated  that  the  optimum  dry  period  was  from  35-40  days 

for  the  Swedish-Friesian  breed,  while  Smith  [1959]  found  it  to 
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be  about  55  days  for  a  constant  calving  interval  of  365  days 
for  10,000  pounds  of  milk  in  American  breeds.   An  examination 
of  the  frequency  distribution  for  previous  days  dry  showed  that 
the  largest  numbers  of  records'  were  for  the  period  from  47  to 
68  days  with  the  largest  number  being  for  59  days.  The  number 
of  days  previously  dry  ranged  from  one  day  to  over  328  days, 
with  the  numbers  increasing  to  59  days  and  decreasing  thereafter. 
[This  discussion  pertains  to  Utah  cows,  as  there  were  more  of 
them  than  those  for  any  other  state.  The  frequency  distributions 
were  considered  on  a  within-state  basis.]   After  studying  the 
frequency  distribution,  it  was  decided  to  subdivide  days 
previous  dry  in  the  manner  indicated  in  Table  7. 

Table  7 
The  data  classes  of  previous  days  dry  and  the  codes  thereof. 
Previous  Days  Dry  Code 

20  days  and  under  1 

21  days  to  40  days  2 
41  days  to  60  days  -3 
61  days  to  80  days  4 
81  days  to  100  days  5 
100  days  and  over  6 


Note  that  previous  days  dry  cannot  be  considered  for  first 
lactations,  and  therefore  they  were  dropped  from  this  part  of 
the  analysis. 

A  study  of  the  frequency  distribution  for  Utah  for  days 
open  showed  that  the  number  of  cows  per  day  was  quite  small 
up  to  69  days  open.   From  69  days  open  to  200  days,  the  number 
of  cows  per  day  is  quite  uniform,  although  it  reaches  a 
maximum  at  92  days  and  is  somewhat  lower  at  the  ends  of  this 
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interval.  This  variable  was  arbitrarily  subdivided  into  the 
intervals  described  in  the  table  below. 

Table  8 
The  data  classes  of  days  open  and  the  codes  thereof. 

Days  Open  Code 

60  days  and  lander  1 

61  days  to  80  days  2 
81  days  to  100  days  3 
101  days  to  120  days  4 
121  days  to  140  days  5 
141  days  to  160  days  6 
161  days  to  180  days  7 
181  days  and  over              8 

The  following  data  classes  for  level  of  milk  production 
were  determined  by  examining  the  frequency  distribution  for  Utah. 

Table  9 

The  data  classes  of  level  of  milk  production  and  the  codes 
thereof. 

Level  of  Milk  Production  Code 

5,999  pounds  and  less  1 

6,000  pounds  to  6,999  pounds         2 

7,000  pounds  to  7,999  pounds        3 

8,000  pounds  to  8,999  pounds  4 

9,000  pounds  to  9,999  pounds         5 

10,000  pounds  to  10,999  pounds  6 

11,000  pounds  to  11,999  pounds  7 

12,000  pounds  to  12,999  pounds  8 

13,000  pounds  to  13,999  pounds  9 

14,000  pounds  to  14,999  pounds  10 

15,000  pounds  to  15,999  pounds  11 

16,000  pounds  to  16,999  pounds  12 

17,000  pounds  to  17,999  pounds  13 

18,000  pounds  to  18,999  pounds  14 

19,000  pounds  to  19,999  pounds  15 

20,000  pounds  and  greater  16 


Level  of  fat  production  was  divided  into  classes  by 
examining  the  mean  and  variance  of  the  frequency  distribution 
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for  Utah.  The  mean  was  found  to  be  approximately  4-60  pounds 
of  fat,  with  a  standard  deviation  of  approximately  100  pounds. 
An  interval  of  roughly  1/3  the  standard  deviation  was  used  and 
the  following  levels  were  determined. 

Table  10 

The  data  classes  of  level  of  fat  production  and  the  codes 
thereof. 

Level  of  Fat  Production  Code 

309  pounds  and  less  1 

310  pounds  to  339  pounds  2 
340  pounds  to  369  pounds  3 
370  pounds  to  399  pounds  4 
400  pounds  to  429  pounds  5 
430  pounds  to  459  pounds  6 
460  pounds  to  489  pounds  V 
490  pounds  to  519  pounds  8 
520  pounds  to  549  pounds  9 
550  pounds  to  579  pounds  10 
580  pounds  to  609  pounds  11 
610  pounds  to  639  pounds  12 
640  pounds  to  669  pounds  13 
670  pounds  to  699  pounds  14 
700  pounds  and  greater  15 

2.2.3.   Classification  of  the  analyses  to  be  performed.  As  can 
be  seen,  if  these  eight  variables  were  analyzed  simultaneously, 
one  would  need  several  million  cells.   For  this  reason,  the 
over-all  analysis  was  broken  into  the  following  analyses. 

1.  A  four-way  analysis  of  Geographic  Area,  Average  Per 
Cent  Fat,  Lactation  Number,  and  Season  for  Geographic  Areas  2 
and  4. 

2.  A  four-way  analysis  of  Geographic  Area,  Average  Per 
Cent  Fat,  Lactation  number,  and  Season  for  Geographic  Areas  3 
and  5- 
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3.  A  two-way  analysis  of  Days  Previous  Dray  and  Days  Open. 

4.  A  one-way  analysis  of  Level  of  Milk  Production. 

5.  A  one-way  analysis  of  Level  of  Pat  Production. 

Two  sets  of  these  analyses  were  performed;  one  for  milk  pro- 
duction, and  one  for  fat  production. 

2.2.4.   Hature  of  the  statistical  problems  involved  and  their 
solutions. 

2.2.4.1.  Before  presenting  the  models  and  the  analyses  for  the 
five  classifications  of  Section  2.2.3,  a  discussion  of  the 
problems  involved  seems  necessary.  To  justify  the  solutions 
obtained  for  some  of  these  problems,  the  discussions  of 
statistical  considerations  of  Section  1.4  will  be  referred  to. 

2.2.4.2.  The  ten  period  approach  vs.  the  305  <3.ay  approach. 
There  are  two  ways  of  expressing  a  lactation  curve.   The  first 
is  to  present  the  ten  monthly  test  days,  or  the  monthly 
production.  The  problem  then  consists  of  studying  ten  equally 
spaced  observations  for  each  cow.  The  second  method  is  more 
closely  associated  with  the  actual  Dairy  Herd  Improvement 
Association  rules  and  practices.  Each  month's  production  is 
based  on  one  test  for  that  month,  and  the  results  are  multiplied 
by  the  number  of  days  in  the  month.  To  make  this  test  meaningful, 
because  of  the  gradual  decline  in  production  over  the  lactation, 

a  centering  date  is  introduced,  causing  the  production  to  be 
estimated  for  15  days  prior  to  that  date  and  12  to  15  days  after 
that  date,  depending  on  the  month  of  test.  The  test  date  should 
fall  within  3  days  of  the  centering  date,  but  in  practice,  the 
difference  is  sometimes  greater.   Because  a  cow  cannot  be 
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analysis,  where  "X  =  days  in  milk"  can  take  on  the  values 
1, — ,  305.   The  first  problem  is  that  there  will  not  be  the 
same  number  of  observations  for  each  record.  The  second  is  that 
the  X  values  will  not  be  the  same  for  all  records.  Despite 
the  fact  that  x,  y  and  N  are  random,  once  they  are 
determined,  X  can  be  considered  a  set  of  fixed  values. 

It  is  difficult  to  use  standard  statistical  analyses  of 
the  8  to  11  observations  in  comparing  the  lactation  estimation 
curves.  Thus,  it  was  decided  to  analyze  differences  between 
the  lactation  estimation  curves  by  using  the  regression  of  Y 
on  X,  where  Y  is  the  ratio  of  the  cumulative  monthly 
production  to  the  total  production  in  305  days.  This  reduces 
the  comparison  of  a  possible  11  observations  per  record  to  one 
observation  per  record.  Fitting  orthogonal  polynomials  to 
factors  representing  the  shape  of  the  lactation  estimation  curve 
indicated  that  when  this  approach  was  used,  it  was  necessary  to 
use  a  fourth  degree  polynomial  to  adequately  represent  the  curve. 
This  would  give  four  coefficients  to  analyze  instead  of  one. 
[Crandall,  19633.  It  would  be  better  to  get  the  curve  in  terms 
of  one  coefficient.   [See  Eao  [1958]  and  Section  1.5. 3]. 
2.2.4.4-.   The  problem  of  the  distribution  of  a  ratio  and  its 
solution.   Since  the  points  on  the  lactation  estimation  curve 
are  defined  to  be 

Y  =  y°tal  305-day  production 
X    Total  production  in  X  days   ' 

and  since  both  the  numerator  and  the  denominator  are  assumed  to 

be  normally  distributed,  it  is  necessary  to  determine  the 
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distribution  of  Y.  From  the  articles  by  Geary  [1930], 
Fieller  [1932],  and  Marsaglia  [1965]  cited  in  Section  1.4.5, 
since  both  the  numerator  and  the  denominator  are  positive  and 
not  close  to  zero,  the  coefficients  of  variation  being  small, 
the  distribution  of  a  function  of  Y  can  be  found  which  is 
approximately  normal.   Because  of  the  difficulty  in  calculating 
the  necessary  standard  deviation  and  correlations,  the  data  were 
not  transformed  to  this  function  which  would  normalize  the 
observations.  Instead  it  was  hoped  that  the  use  of  the  trans- 
formation of  Section  2.2.4.5  scud,   the  use  of  means  of  observations 
[Section  2.2.4.6]  would  "smooth"  the  data  and  cause  them  to 
approach  a  normal  distribution. 

2.2.4.5.   The  problem  of  a  transformation  and  its  solution. 
The  lactation  estimation  curve  is  distinctly  curvilinear,  and 
thus  if  linear  regression,  as  discussed  above  [Section  2.2.4.3] 
were  used,  a  poor  fit  would  result.   Instead  the  curvilinear 
model  was  assumed  to  be 

Y  =  A[305/X]C.  [5] 

In  order  to  fit  this  equation,  the  log  of  both  sides  was  taken 
to  give 

log10  Y  =  log10  A  +  C  log10  [305/X].  [6] 

This  equation  was  fitted  for  each  record  using  standard 
regression  techniques  [soe  Appendix  I.]   The  analysis  of 
variance  for  regression  on  a  small  portion  of  the  data 
indicated  extremely  high  significance.  When  two  cases  of 
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actual  data  were  plotted  on  log-log  paper,  [Figure  1]  there 
was  still  a  slight  curvelinear  effect  after  the  transformation. 
However,  when  overall  factors  obtained  by  this  method  were 
compared  with  factors  in  current  use,  it  was  found  that  there 
was  a  very  close  fit. 

Figure  1  illustrates  that  heterogeneity  among  a  group  of 
fitted  regression  lines  will  cause  the  largest  differences  to 
occur  early  in  the  lactation,  and  the  lines  to  converge  to  one 
at  305  days.   [Note  the  difference  between  Equation  (5)  and  that 
of  Figure  1.  The  regression  coefficients  were  all  calculated 
using  the  equation  in  Figure  1,  which  yielded  a  negative  co- 
efficient for  all  records.  By  omitting  the  minus  sign,  for  all 
coefficients,  one  has  the  same  results  as  would  have  been  obtained 
had  Equation  (5)  been  used.]  The  reason  why  this  heterogeneity 
exists  is  because  Yj  does  not  have  uniform  variance  for  all 
values  of  X.  When  X  is  small,  the  variance  of  Yx  is 
extremely  large,  and  as  X  approaches  305  days,  the  variance 
of  Yx  approaches  zero.  This  effect  of  the  variance  violates 
one  of  the  assumptions  of  covariance  analysis,  i.e.,  that  of 
homogeneous  variance  along  the  regression  line,  but  does  not 
affect  the  normality  of  the  regression  coefficients.   Thus 
analyses  of  variance  can  be  performed  using  the  regression 
coefficients  as  observations,  but  analysis  of  covariance 
techniques  cannot  be  used. 

When  the  regression  coefficients  were  examined,  it  was 
found  that  the  values  of  A  [Equation  (5)]  were  practically  one. 
thus  this  term  was  not  considered  in  later  studies.  Actually, 
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a  more  accurate  model  would  have  been  obtained  by  ignoring  the 
mean  and  recalculating  the  regression  coefficients,  thus  giving 

Y  =  [305AJC.  C7] 

In  this  case,   C  =  E  X'Y'/E  X'2,  where  X'  =  log  [3O5/X]   and 

p  

Y'  =  log  Y,   instead  of  C  =  E  x'y'/E  x'  ,   where  x'  ■  [X1  -x'J 

and  y  ■  [Y1  -y'].  While  this  difference  is  unimportant  with 

respect  to  the  regression  coefficient,  C,   it  does  make  a 

difference  if  covariance  analysis  is  to  be  used.  The  fact  that 

one  must  subtract  [E  X']  /N  from  the  denominator  of  the 

standard  error  makes  a  larger  standard  error  when  the  mean  is 

considered,  incorrectly,  in  the  model,  than  if  the  mean  were  not 

considered.  For  these  reasons,  it  was  decided  to  perform 

analysis  of  variance  techniques  with  the  regression  coefficients 

for  each  record  as  the  observations,  rather  than  applying 

covariance  techniques. 

2.2.4.6.  The  problem  of  unequal  subclass  frequencies  and  its 

solution.   Because  of  unequal  subclasses  involved  in  this  study, 

least  squares  techniques  seemed  applicable,  however,  none  of 

the  least  squares  computer  programs  available  at  either  Brigham 

Young  University  [IBM  7040]  or  the  University  of  Washington 

[IBM  7094]  would  handle  as  large  an  analysis  as  the  first  two 

described  in  Section  2.2.3.   ^n  looking  for  an  alternate 

method,  it  was  decided  to  use  the  method  of  unweighted  means 

[Yates,  1934].   A  second  consideration  for  this  particular 

analysis  was  that  it  was  felt  that  the  means  of  the  regression 

coefficients  would  be  more  nearly  normal.  For  this  method, 

the  means  were  analyzed  according  to  a  standard  orthogonal 
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factorial  design  for  the  first  three  analyses  described  in 
Section  2.2.3,  and  a  one-way  analysis  of  variance  with  unequal 
subclass  frequencies  was  used  for  the  last  two  analyses  described. 
In  the  factorial  analyses,  the  highest  order  interaction  was 
used  as  an  estimate  for  error. 

2.2.5.   The  linear  statistical  models  for  the  five  analyses. 
All  effects  are  considered  fixed,  with  the  exception  of  the 
error  term. 

For  the  first  analysis  for  areas  2  and  4-,  the  regression 
coefficients  for  milk  and  fat  were  each  assumed  to  be  described 
by  the  following  linear  model • 

^iaki  "  v-  +  ai  +  b 3   +  ck  +  di  +  abia  +  acik  +  adn  +  bcjk  + 

♦  bddl  ♦  cdkl  *  abcijk  +  abdi;jl  ♦  acdikl  +  bcdjkl  +  e.  .^ 

i  =  2,4;      k  -  1,...,  6; 

i   =  2,...,  4-;  1  =   1,...,  4;   for  all  yiikl  in  areas  2  and  4-. 
where 
yijkl  =  tixe  TeS^essi-°n   coefficient  in  the  ith  geographic  area, 

in  the  j   per  cent  fat  level,  in  the  kt]l  lactation,  and 
in  the  1   season 
|j.     =  the  overall  regression  coefficient  for  the  population, 
when  equal  frequencies  exist  in  all  subclasses. 
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=  the  effect  of  the  i   geographic  area. 


b  .    =  the  effect  of  the  oth  per  cent  fat  level. 


c,     =  the  effect  of  the  kt]l  lactation. 


d]_    =  the  effect  of  the  1 

abin  =  tile  e££ea*  o£   tlle  ij   subclass  of  geographic  area  and 


32 


per  cent  fat  after  the  average  effects  of  the  i 
geographic  area  and  the  j   per  cent  fat  level  have 
been  removed. 

=  the  effect  of  the  ik   subclass  of  geographic  area  and 
lactation  number  after  the  average  effects  of  the  i 
geographic  area  and  the  k   lactation  have  been  removed. 

=  the  effect  of  the  il   subclass  of  geographic  area  and 
season  after  the  average  effects  of  the  i   geographic 


be..   =  the  effect  of  the  jk   subclass  of  per  cent  fat  and 

lactation  after  the  average  effects  of  the  j   per  cent 
fat  level  and  the  k   lactation  have  been  removed. 

bd.,   =  the  effect  of  the  jl   subclass  of  percent  fat  and 


Ml 


cdkl 


season  after  the  average  effects  of  the  j   per  cent 


after  the  average  effects  of  the  k  ,  lactation  and  the 


abo  ..  ■  the  effect  of  the  ijk   subclass  of  geographic  area, 
per  cent  fat,  and  lactation  number  after  the  average 
effects  of  the  i   geographic  area,  the  j   per  cent 
fat  level,  the  .k   lactation,  the  ij   subclass  of 
geographic  area  and  per  cent  fat,  the  ik   subclass  of 
geographic  area  and  lactation  number,  and  the  jk 
subclass  of  per  cent  fat  and  lactation  number  have  been 
removed. 
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abd. .,  =  the  effect  of  the  ijl   subclass  of  geographic  area, 
per  cent  fat,  and  season  after  the  average  effects  of 
the  i   geographic  area,  the  j   per  cent  fat  level, 
the  1   season,  the  ij   subclass  of  geographic  area 
and  per  cent  fat,  the  il   subclass  of  geographic  area 
and  season,  and  the  jl   subclass  of  per  cent  fat  and 
season  have  been  removed. 

acd..  ,  =  the  effect  of  the  ikl   subclass  of  geographic  area, 

lactation  number,  and  season  after  the  average  effects 
of  the  i   geographic  area,  the  k   lactation  number, 
the  1   season,  the  ik   subclass  of  geographic  area 
and  lactation  number,  the  ik   subclass  of  geographic 
area  and  season,  and  the  kl   subclass  of  lactation 
number  and  season  have  been  removed. 

bcd.kl  =  the  effect  of  the  dkl   subclass  of  per  cent  fat, 

lactation  number,  and  season  after  the  average  effects 
of  the  o   per  cent  fat  level,  the  k   lactation,  the 
1   season,  the  jk   subclass  of  per  cent  fat  level 
and  lactation  number,  the  jl   subclass  of  per  cent 
fat  level  and  season,  and  the  kl   subclass  of 
lactation  number  and  season  are  removed. 

ei-jjcl  =  the  residual  error,  the  amount  of  variation  not 

accounted  for  by  the  above  effects, [the  highest  order 

interaction  assumed  to  be  zero] .  These  errors  are 

p 
assumed  to  be  NID[0,  cr     ] . 

The  various  effects  are  expressed  as  deviations  from  the  mean, 
and  therefore,   £  a.  ■  2  b.  =  S  ck  =  £  d,  =  0,   and 
2  abj.  -  2  ab.j.  =  ...  =  E  cdkl  =  E  cdkl  =  0.  Similar 

restrictions  are  placed  on  the  three-way  interactions.  Because 
all  the  effects  in  this  study  are  assumed  fixed,  estimation  of 
the  constants  for  the  levels  of  the  significant  effects  were 
made  to  be  used  in  estimating  total  production  from  part 
lactation  records. 
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For  the  second  analysis,  for  areas  3  and  5,  the  model  is 
exactly  the  same,  except  that  the  subscripts,   i,...,  1.  are 
over  different  values,   i  =  3,  5;   0  =  2,...,  6;  k  =  1,...,  6; 


and 


1.1,..,,  4|   for  all  yidkl  in  areas  3  and  5- 
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above  assumptions  also  hold  in  this  model. 

For  the  third  analysis,  for  days  dry  and  days  open,  the 
regression  coefficients  for  milk  and  fat  were  each  assumed  to 
be  described  by  the  following  model 

yi;jk  - 1» +  ai +  b.j +  eio 

i  =  1,...,  6  and  j    =   1,...,  8 
where 
y.   ■  the  regression  coefficient  in  the  i   class  of  days 

previous  dry  and  the  j   class  of  days  open. 
H   =  the  overall  regression  coefficient  for  the  population  when 

equal  frequencies  exist  in  all  subclasses. 

a.  =  the  effect  of  the  i   class  of  days  previous  dry. 

b.  =  the  effect  of  the  oth  class  of  days  open. 

u 

e. .  =  random  errors  not  accounted  for  by  the  other  terms  in 

the  model. 

Since  the  interaction  is  assumed  to  be  zero,  the  error 

term  is  represented  by  the  two-way  interaction,  and  the  errors 

o 
expressed  thereby  are  assumed  NID  [0,  a~e   ] .   It  should  be  noted 

that  the  y- j  are  from  the  entire  population  of  areas  2  and  4, 

as  well  as  3  and  5,  excluding  the  y^  for  which  days  carried 

calf  equal  zero  and  those  for  which  the  record  is  a  first 

lactation. 
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For  the  fourth  analysis,  for  level  of  milk  production,  the 
regression  coefficients  for  milk  and  fat  were  assumed  to  he 
expressed  hy  the  following  linear  model 

y±d  =  i» +  ai +  eio 

i  =  1 , . . . ,  16  and  j  =  1 , . . . ,  tu 
where 

y .  .  =  the  i       regression  coefficient  in  the  i   level. 
|i   =  the  overall  regression  coefficient  for  the  population 
when  equal  frequencies  exist  in  all  subclasses. 

&,      =  the  effect  of  the  i   level  of  milk  production. 

p 
e.  .  =  random  errors,  which  are  NID[0,  <r     ]. 

Since  the  a.   are  expressed  as  deviations  about  the  mean, 

S  a^  =  0.   It  should  be  noted  that  the  y^.,  are  here  from  the 

entire  population,  areas  2,  3,  4-,  and  5. 

For  the  fifth  analysis,  for  level  of  fat  production,  the 
regression  coefficients  for  milk  and  fat  were  assumed  to  be 
described  by  the  following  model 
y.  .  =  u  +  a.  +  e .  . 

i  ■  1,...,  15  and  j  =  1,...,  n^ 
where 

yi1  =  the  ^   regression  coefficient  in  the  i   level. 
\l       =  the  overall  regression  coefficient  for  the  population  when 

equal  frequencies  exist  in  all  subclasses. 
ai  =  the  effect  of  the  i   level  of  fat  production, 
e.^j  =  random  errors,  which  are  NID[0,  «r -J. 

The  same  comments  hold  here  as  for  the  fourth  analysis  above. 


CHAPTER  3 
RESULTS  AND  DISCUSSION 

3.1-  Results  of  the  analyses  of  variance  and  discussion.  The 
results  of  analyses  1  through  5  for  milk  and  fat  will  be 
presented  in  the  following  tables. 

Table  11 

Analysis  of  regression  coefficients  for  factors  for  estimating 
total  milk  production  from  cumulative  milk  production  for 
geographic  areas  2  and  4. 


Source 

df 

Sums  of  Squares 

Mean  Squares 

F 

Geographic  Area 

1 

0.00000 

0.00000 

0.00000 

Per  Cent  Fat 

2 

0.01489 

O.OO745 

17.32558** 

Lactation  Number 

5 

0.09036 

0.01807 

42.02326** 

Season 

3 

0.03269 

0.01090 

25.34884** 

GA  x  PF 

2 

0.00051 

0.00026 

0.60465 

GA  x  LN 

5 

0.00561 

0.00112 

2.60465* 

GA  x  S 

3 

0.00098 

O.OOO33 

0.76744 

PF  x  ln 

10 

0.00388 

0.00039 

0.90698 

PF  x  S 

6 

0.01487 

0.00248 

5.76744-** 

LN  x  s 

15 

0.01056 

0.00070 

1.62791 

GA  x  pp  x  LN 

10 

0.00645 

0.00064 

1.48837 

GA  x  PF  x  s 

6 

0.00174 

.0. 00029 

0.67442 

GA  x  LN  x  S 

15 

0.00765 

0.00051 

1.18605 

PF  x  LN  x  S 

30 

0.01984 

0.00066 

1.53488 

Residual 

30 

0.01292 

0.00043 

Total 


143 


0.22296 


Note  that  wherever  [*] 
"significant  at  the  a 


is  used  in  an  F-column,  it  means 
=  O.05  level",  and  wherever  [**]  is 

used  in  an  F-column,  it  means  "significant  at  the  a  =  0.01 

level. 
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Table  12 

Analysis  of  regre 

ssion 

coefficients  for 

factors  for  e: 

stimating 

total  milk  production 

from  cumulative  milk  production 

for 

geographic  areas 

5  and 

5- 

Source 

df  i 

Sums  of  Squares 

Mean  Squares 

F 

Geographic  Area 

1 

0.00113 

0.00113 

3.53125 

Per  Cent  Fat 

4 

0.05572 

0.01393 

43.53125** 

Lactation  Number 

5 

0.11283 

0.02257 

70.53125** 

Season 

3 

0.07059 

0.02353 

73.53125** 

GA  x  PF 

4 

0.00412 

0.00103 

3. 21875* 

GA  x  LN 

5 

0.00228 

0.00046 

1.43750 

GA  x  S 

3 

0.00021 

0.00007 

0.21875 

PF  x  LN 

20 

O.OI383 

0.00069 

2.15625* 

PF  x  S 

12 

0.01114 

0.00093 

2.90625** 

LN  x  s 

15 

0.00160 

0.00011 

0.3*375 

SA  x  PF  x  LN 

20 

0.00499 

0.00025 

0.78125 

SA  x  pp  x  s 

12 

0.00285 

0.00024 

0.75000 

SA  x  LN  x  S 

15 

0.00853 

0.00057 

1.78125 

PF  x  LN  x  S 

60 

0.01413 

0.00024 

0.75000 

Residual 

60 

0.01916 

0.00032 

Total 

239 

O.323II 

Table  13 

p 

Analysis  of  regre 

ssion 

coefficients  for 

factors  for  e: 

stimating 

total  milk  production 

from  cumulative  milk  production 

for  all 

areas,  excluding 

first 

lactation  records 

1  and  those  for  which 

there  are  no  days 

carried  calf. 

Source 

df 

Sums  of  Squares 

Mean  Squares 

F 

Days  Previous  Dry 

5 

0.18584 

0.03717 

1.31204 

Days  Open 

7 

0.18667 

0.02667 

0.9*140 

Residual 

35 

0. 99163 

0.02833 

Total 

47 

1.36415 
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Table  14 

Analysis  of  regression  coefficients  for  factors  for  estimating 
total  milk  production  from  cumulative  milk  production  for  all 
records  in  all  areas. 


Source 

df 

Sums  of  Squares 

Mean  Squares    P 

level  of 
Prod. 

Residual 

Milk 

15 
18,529 

1.95564 

118.46874 

0.13038    20.40378' 
0.00639 

Total 

18,544* 

120.42438 

*Note  that  there  are  18,541  records  used  in  the  study,  while 
18,545  records  were  used  in  this  analysis.  This  is  because  the 
sort  routine  used  to  set  up  this  analysis  filled  in  the  four 
remaining  positions  in  the  last  block  of  tape  with  four  valid 
records  from  the  last  cell.  Since  there  were  175  records  in 
the  last  cell,  it  was  felt  that  the  effects  of  this  duplication 
would  be  negligible. 


Table  15 

Analysis  of  regression  coefficients  for  factors  for  estimating 
total  milk  production  from  cumulative  milk  production  for  all 
records  in  all  areas. 


Source 


df   Sums  of  Squares   Mean  Squares 


Level  of  Pat 
Prod. 


Residual 


14 
18,530 


0.52779 
119.14661 


O.O377O 
0.00643 


5.86314** 


Total 


18,544*   119.67440 


*Note  that  the  discussion  following  Table  14  applies  here 
also,  since  there  were  227  records  in  the  last  cell  for  this 
analysis. 
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Table  16 

Analysis  of  regression 

coefficients  for 

factors  for 

estimating 

total  fat  production  from  cumulative  fat 

;  production 

for 

geographic  areas 

2  and 

4. 

Source 

df 

Sums  of  Squares 

Mean  Squares    F 

Geographic  Area 

1 

0.00263 

0.00263 

5.4-9717* 

Per  Cent  Fat 

2 

0.00607 

0.00304 

6.33333** 

Lactation  Number 

5 

0.13302 

0.02660 

55.^1667** 

Season 

3 

0.00663 

0.00221 

4.60417** 

GA  x  pf 

2 

0.00286 

0.00143 

2.97917 

GA  x  LN 

5 

0.00418 

0.00084 

1.75000 

GA  x  s 

3 

0.00212 

0.00071 

1.47917 

PF  x  LN 

10 

0.01006 

0.00101 

2.10417 

PF  x  s 

6 

0.01252 

0.00209 

4.35417** 

LN  x  s 

15 

0.00812 

0.00054 

1.12500 

GA  x  pf  x  ln 

10 

0.00810 

0.00081 

I.68750 

GA  x  PF  x  s 

6 

0.00402 

0.00067 

1.39853 

GA  x  LN  x  s 

15 

O.OO783 

0.00052 

1.08333 

PF  x  LN  x  S 

30 

0.02747 

0.00092 

1. 91666* 

Residual 

30 

0.01455 

0.00048 

Total 

143 

0.25018 

Table  17 

Analysis  of  regression 

coefficients  for 

factors  for 

estimating 

total  fat  production  from  cumulative  fat 

production 

for 

geographic  areas 

3  and 

5- 

Source 

df 

Sums  of  Squares 

Mean  Squares    F 

Geographic  Area 

1 

0.00094 

0.00094 

1.84314 

Per  Cent  Fat 

4 

0.00874 

0.00219 

4.29412** 

Lactation  Number 

5 

0.20104 

0.04021 

78.84314** 

Season 

3 

0.00476 

0.00159 

3.H765* 

GA  x  pf 

4 

0.01281 

0.00320 

6.27451** 

GA  x  LN 

5 

0.00206 

0.00041 

0.80392 

GA  x  s 

3 

0.00031 

0.00010 

0.19608 

PF  x  LN 

20 

0.00883 

0.00044 

0.86275 

PF  x  S 

12 

0.00457 

O.OOO38 

O.74510 

LN  x  S 

15 

0.0061 3 

0.00041 

0.80 392 

GA  x  pf  x  LN 

20 

O.OO717 

0.00036 

O.70588 

GA  x  pf  x  S 

12 

0.00363 

0.00030 

0.58824 

GA  x  LN  x  S 

15 

0.00951 

0.00063 

1.23529 

PF  x  LN  x  S 

60 

0.01801 

0.00030 

0.58824 

Residual 

60 

0.03042 

0.00051 

Total 

239 

0. 31893 

"  ■ 
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Table  18 

Analysis  of  regression  coefficients  for  factors  for  estimating 
total  fat  production  from  cumulative  fat  production  for  all 
areas,  excluding  first  lactation  records  and  those  for  which 
there  are  no  days  carried  calf. 


Source 


df   Sums  of  Squares   Mean  Squares 


Days  Previous  Dry  5  0.14375 
Days  Open  7  0.18927 
Residual         35     0.99977 


0.02875 
0.02704 
0.2856 


1.00067 
0.94678 


Total 


47 


1.33279 


Table  19 

Analysis  of  regression  coefficients  for  factors  for  estimating 
total  fat  production  from  cumulative  fat  production  for  all 
records  in  all  areas. 


Source 


df   Sums  of  Squares   Mean.  Squares 


Level  of  Milk 

Prod.        15 

1.63407 

Hesidual      18,529 

I47.3078O 

0.10894 
0.00795 


13.70314** 


Total         18,544*   148.94187 
*  See  the  note  for  Table  14. 


Table  20 

Analysis  of  regression  coefficients  for  factors  for  estimating 
total  fat  production  from  cumulative  fat  production  for  all 
records  in  all  areas. 


Source 


df   Sums  of  Squares   Mean  Squares    F 


Level  of  Fat 
Prod. 

Residual 


14 
18,530 


1.05059 
147. 86122 


0.07504 
0.00798 


9.40351" 


Total         18,544*   148.91181 
*  See  the  note  for  Table  15. 


0.00798 
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A  study  of  the  above  tables  indicates  that  Per  Cent  Pat, 
Lactation  Number,  and  Season  are  highly  significant  for  all 
areas  for  both  milk  and  fat  factors,  with  the  exception  of 
Season  for  the  analysis  of  fat  for  areas  3  and  5.  For  the 
latter  analysis  the  effect  of  Season  is  still  significant. 
Thus  Breed,  Age,  and  Season  all  affect  the  shape  of  the  lactation 
estimation  curve.   Geographic  Area  was  found  to  be  significant 
only  in  the  case  of  the  analysis  for  fat  in  areas  2  and  4. 
The  marginal  means  for  Geographic  Area  in  this  case  of  milk  are 
0.87916  and  0.87897,  for  areas  2  and  4,  respectively.  As  can 
be  seen,  these  are  practically  the  same  and  hence  the  zero  sum 
of  squares  in  Table  11.  Prior  to  combining  areas  1  and  2,  an 
analysis  was  performed  using  areas  1,  2,  and  4.   In  this  analysis, 
the  effect  of  Geographic  Area  was  found  to  be  highly  significant 
for  milk,  but  not  significant  for  fat.   In  the  analysis  for  fat 
in  this  case,  the  effect  of  Per  Cent  Fat  was  also  found  to  be 
non-significant.   In  spite  of  this,  areas  1»  and  2  were  combined 
because  there  were  so  few  records  from  these  areas.   It  should 
be  remembered  that  the  analyses  for  Tables  11,  12,  16,  and  17 
[as  well  as  13  and  18]  are  for  the  cell  means  only,  thus,  the 
means  for  the  cells  in  areas  1  and  2  are  considered  in  the 
analysis  to  have  the  same  weight  as  those  for  area  4.   Referring  - 
to  Table  1,  it  can  be  seen  that  there  are  approximately  4  times 
as  many  records  in  area  2  [second  three  states]  as  in  area  1 
[first  three  states]  and  5  times  as  many  records  in  area  4 
[New  Mexico-Arizona]  as  in  area  1.  When  the  first  two  areas 
are  combined,  the  numbers  of  records  in  the  two  groups  to  be 
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compared  are  approximately  equal.  Because  of  this  disparity 
in  weighting,  then,  it  was  decided  to  use  the  combined  analysis 
given  in  Tables  11  and  16.   It  might  he  noted,  however,  that  the 
regression  coefficients  in  area  1  were  consistently  higher  than 
those  for  area  2. 

There  is  a  highly  significant  interaction  between  Per 
Cent  Fat  and  Season  for  milk  and  fat  in  areas  2  and  4  and 
for  milk  in  areas  3  and-  5»  which  seems  to  indicate  that  different 
breeds  [or  cows  producing  different  fat  percentages]  react 
differently  to  different  seasonal  conditions.  This  effect  also 
appears  in  both  milk  and  fat  analyses  for  areas  3  an<i  5  with 
respect  to  the  interaction  between  Breed  [per  cent  fat]  and 
Geographic  Area.   In  this  particular  case,  the  significant 
interaction  indicates  that  cows  with  different  percentages  of 
fat  will  have  a  different  lactation  estimation  curve  in 
Colorado  than  those  in  Utah  or  Nevada.  For  fat,  this  Area  by 
Breed  interaction  effect  seems  actually  gre.ater  than  the  effect 
of  Season  or  Breed  [see  Section  3.2]  on  the  lactation  curve. 
One  might  also  note  that  in  areas  2  and  4  for  milk,  there  is  a 
significant  interaction  between  Geographic  Area  and  Lactation 
Number.  Also  in  areas  3  and  5  for  milk,  there  is  a  significant 
interaction  between  Lactation  Number  and  Per  Cent  Fat.   In  the 
analysis  for  fat  in  areas  2  and  4,  there  is  a  significant  three- 
way  interaction  between  Per  Cent  Fat,  Lactation  Number  and 
Season.   The  effects  of  these  interactions  will  be  looked  at 
more  closely  in  the  section  dealing  with  the  analysis  of  the 
means  [see  Section  3.3]. 
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An  examination  of  Tables  13  and  18  indicates  that  neither 
Days  Previous  Dry  nor  Days  Open  seem  to  have  any  significant 
effect  on  the  shape  of  the  lactation  estimation  curve  for 
either  milk  or  fat. 

Tables  14,  15,  19  and  20  indicate  that  both  Level  of 
Milk  Production  and  Level  of  Pat  Production  are  highly  signifi- 
cant for  both  milk  and  fat.  The  fact  that  there  is  an  extremely 
large  number  of  degrees  of  freedom  associated  with  the  within 
or  residual  sum  of  squares  could  indicate  that  actually  there 
is  a  very  small  difference  with  a  very  sensitive  test.  Since 
these  tests  of  significance  converge  to  significance  as  the 
degrees  of  freedom  become  infinite,  this  could  very  well  be 
the  case  here.  The  means  of  these  analyses  will  be  looked  at 
in  greater  detail  in  Section  3.3. 

3-2.  Results  of  the  analyses  of  the  components  of  variance  and 
discussion.  Although  this  study  is  concerned  with  the  fixed 
effects  model  and  thus,  one  is  not  primarily  concerned  with  the 
estimation  of  the  variance  components,  it  was  felt  that  much 
valuable  information  could  be  obtained  by  studying  the  "variances" 
of  the  constants  involved  in  this  finite  population  to  determine 
the  relative  importance  of  the  significant  effects,  i.e.,  the 
proportion  of  the  total  variance  for  which  they  account.   Since 
a  fixed  effects  model  is  being  assumed,  the  expected  mean 
squares  were  calculated,  equated  to  their  respective  mean 
squares,  and  solutions  obtained  for  estimates  of  the  components 
of  variance.  The  entire  analysis  is  based  on  the  means  of  the 
cells,  hence,  the  analyses  of  the  components  of  variance  are 
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considered 

as  the 

orthogonal  case  with  one 

observation  per 

cell. 

Belovv 

are  presented  the  four  four-way 

analyses  in  terms 

of  their  c 

omponents  of 

variance. 

Table  21 

• 
analysis  of  milk 

Analysis  of  components 

of  variance  for  the 

coefficients  for 

geographic  areas  2  and  4. 

[See  Table  11.] 

Source 

df 

Component  of  Variance   Coded     % 

Geographic 

Area 

1 

0.0000000 

0     0.000 

Per  Cent  Fat 

2 

0.0001462 

1462**   7.065 

Lactation 

Number 

5 

0.0007350 

7350**   35- 516 

Season 

3 

0.0002908 

2908**   14.052 

GA  x  PF 

2 

0.0000000 

0     0.000 

GA  x  M 

5 

0.0000575 

575*    2.778 

GA  x  s 

3 

0.0000000 

0     0.000 

PF  x  LN 

10 

0.0000000 

0     0.000 

PF  x  S 

6 

0.0001708 

1708**   8.253 

LN  x  s 

15 

0.0000450 

450     2.174 

GA  x  PF  x 

LN 

10 

0.0000525 

525     2.537 

GA  x  pp  x 

S 

6 

0.0000000 

0     0.000 

GA  x  LN  x 

S 

15 

0.0000267 

267     1.290 

PF  x  LN  x 

S 

30 

0.0001150 

1150     5-557 

Residual 

30 

0.0004300 

4300    20.778 

Total 

14J 

0.0020695 

20695   100.000 

Note  that 

negative  components  of  variance 

are  considered  to  be 

zero  and  are  reproduced  here  as  such. 
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Table  22 

Analysis  of  - 

components 

of  variance 

for  the 

analysis  1 

of  milk 

coefficients 

for 

geographic  areas  3 

and  5» 

[See  Table  12]. 

Source 

df 

Component  of 

Variance   Coded 

% 

Geographic  Area 

1 

0.0000068 

68 

0.391 

Per  Cent  Fat 

4 

0.0002835 

2835** 

16. 320 

Lactation  Number 

5 

0.0005563 

5563** 

32.025 

Season 

3 

0.0003868 

3868** 

22.267 

GA  x  PF 

4 

0.0000296 

296* 

1.704 

GA  x  LN 

5 

0.0000070 

70 

0.404 

GA  x  S 

3 

0.0000000 

0 

0.000 

PF  x  LN 

20 

0.0000463 

463* 

2.665 

PF  x  s 

12 

0.0000508 

508** 

2.924 

LN  x  S 

15 

0.0000000 

0 

0.000 

GA  x  pp  x  LN 

20 

0.0000000 

0 

0.000 

GA  x  PF  x  S 

12 

0.0000000 

0 

0.000 

GA  x  LN  x  s 

15 

0.0000500 

500 

2.878 

PF  x  LN  x  S 

60 

0.0000000 

0 

0.000 

Residual 

60 

0.0003200 

3200 

18.422 

Total 

239 

0.0017371 

17371 

100.000 

See  note  on  Table  21. 

Table  23 

Analysis  of  components 

of  variance  for  the 

analysis  of  fat 

coefficients 

for 

geographic  areas  2 

and  4. 

[See  Table  16]. 

Source 

df 

Component  of 

Variance  Coded 

% 

Geographic  Area 

1 

0.0000299 

299* 

1.278 

Per  Cent  Fat 

2 

O.OOOO533 

533** 

2.278 

Lactation  Number 

5 

0.0010883 

IO883** 

46.513 

Season 

3 

0.0000481 

481** 

2.056 

GA  x  pf 

2 

0.0000395 

395 

1.688 

GA  x  LN 

5 

0.0000300 

300 

1.282 

GA  x  s 

3 

0.0000128 

128 

0.547 

PF  x  LN 

10 

0.0000662 

662 

2.829 

PF  x  s 

6 

0.0001342 

1342** 

5-736 

LN  x  S 

15 

0.0000100 

100 

0.427 

GA  x  pf  x  LN 

10 

0.0000825 

825 

3-526 

GA  x  pf  x  s 

6 

O.OOOO3I7 

317 

1-355 

GA  x  LN  x  s 

15 

O.OOOOI33 

133 

0.567 

PF  x  LN  x  S 

30 

0.0002200 

2200* 

9.403 

Residual 

30 

0.0004800 

4800 

20.515 

Total 

14-3 

0.0023398 

23398 

100.000 

46 


Table  24 

Analysis  of  components  of  variance  for  the  analysis  of  fat 
coefficients  for  geographic  areas  3  and  5.   [See  Table  17] • 


Source 


df   Component  of  Variance   Coded 


Geographic  Area 

1 

0.0000036 

36 

0.034 

Per  Cent  Fat 

4 

O.OOOO35O 

350** 

0.329 

Lactation  Number 

5 

0.0099250 

99250** 

93.388 

Season 

3 

0.0000180 

180* 

0.169 

GA  x  py 

4 

0.0001121 

1121** 

1.055 

GA  x  LN 

5 

0.0000000 

0 

0.000 

GA  x  s 

3 

0.0000000 

0 

0.000 

PF  x  LN 

20 

0.0000000 

0 

0.000 

PF  x  S 

12 

0.0000000 

0 

0.000 

LN  x  S 

15 

0.0000000 

0 

0.000 

GA  x  pf  x  LN 

20 

0.0000000 

0 

0.000 

GA  x  p?  x  S 

12 

0.0000000 

0 

0.000 

GA  x  LN  x  s 

15 

0.0000240 

240 

0.226 

PF  x  LN  x  S 

SO 

0.0000000 

0 

0.000 

Residual 

60 

0.0005100 

5100 

4.799 

Total  239     0.0106277 

See  note  on  Table  21. 


106277   100.00 


From  these  tables,  it  can  be  seen  that  the  component  of 
variance  for  Lactation  Number  accounts  for  more  of  the  variation 
than  any  other  source  for  the  entire  set  of  analyses  for  fat. 
In  areas  3  and  5,  it  accounts  for  over  93%  of  the  variation  for 
fat.  Lactation  Number  accounts  for  a  much  larger  proportion  of 
the  variation  in  fat  than  it  does  for  milk.  There  is  a  larger 
difference  between  proportions  for  milk  and  fat  in  areas  3  and  5 
than  in  areas  2  and  4.  Per  Cent  Fat  and  Season  account  for  a 
larger  portion  of  the  variation  for  milk  than  they  do  for  fat. 
It  would  seem,  therefore,  that  Per  Cent  Fat  and  Season  have  a 
greater  effect  on  milk  coefficients  than  they  do  on  fat 
coefficients.   In  other  words,  Per  Cent  Fat  [breed]  and  Season 
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exert  a  larger  effect  on  the  lactation  estimation  curves  for 

milk  than  they  do  on  those  for  fat,  and  Lactation  Number  exerts 

a  larger  effect  on  the  Lactation  estimation  curve  for  fat  than 

it  does  for  milk,  although  the  over-all  affect  of  Lactation 

Number  is  larger  than  for  any  other  variables  for  both  milk 

and  fat. 

3.3.   Results  of  the  analysis  of  means  and  discussion.   Graphs 

of  the  significant  variables  and  interactions  are  presented  below. 

Tables  of  comparisons  of  means  using  Tukey's^,  hsd  test  for  these 

effects  are  given  below. 

Table  25 

Table  of  means  for  the  significant  effects  in  Table  12, 
geographic  areas  2  and  4  for  milk,  using  Tukey's  hsd  test. 


hsd  =  0.01045. 


4  0.86468 


Per  Cent  Fat 

0.00299, 

Qq.05  =  5.«, 

2 
0.88639 
0.02171* 

3 
0.88612 
0.02144* 

0.00027 


If 

The  reason  for  the  use  of  Tukey's  test  instead  on  one  of  the 
other  available  tests  for  multiple  comparisons  of  means  is 
explained  in  Section  1.4. 3.   It  was  used  because  of  the  equal 
subclasses  involved  in  the  analyses.   Had  an  analysis  of  the 
means  for  Level  of  Production  for  milk  or  fat  been  given, 
Sheffe's  test  would  have  been  the  ideal  test  to  use,  because 
of  the  unequal  subclass  numbers  involved.   It  should  be  noted 
that  the  a-level  [Type  I  error  probability]  here  is  0.05  for 
each  level  and  not  for  the  over-all  effect  considered.   Thus 
in  Table  26,  the  probability  of  an  error  of  Type  I  is  0.05  for 
Per  Cent  Fat  for  the  difference  between  level  4  and  level  2. 
It  is  also  0.05  for  each  of  the  other  differences  represented 
for  the  levels  of  Per  Cent  Fat.   If  all  three  comparisons 
for  Per  Cent  Fat  are  taken  together,  the  alpha-level  is 
considerably  greater. 
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Table   25 

(Continued) 

Table   of  means   for  the   significant   effects   in  Table   12, 
geographic   areas   2  and  4  for  milk,    using  Tukey's  hsd  test. 

Lactation  Number 

ss  =   0.00423,      q0> 

05  "  4.30,      hsd  «  0. 01819. 

3 
4 
2 
5 
6 

0.86261 
0.86292 
0.86968 
0.87085 
0.87402 

1 
0.93431* 
O.O7T7G* 
O.O7I39* 
0.06463* 
0.06346* 
0.06029* 

6                    5                    2                    4 
0.87402        0.87085        0.86268        0.86292 
U70TI4T       0T0~08~24       O'l'O'OVO?       O.OOO3I 
0.01110       O.OO793       0.00676 
0.00434       0.00117 
O.OO3I7 

Season 

ss   =   0.00346,      qQ> 

05   =   3.84,      hsd  =  0.01329. 

3 

2 
4 

0.86305 
0.86708 
0.88546 

1 
0. 90067 
0.03765* 
0.03359* 
0.01521 

4                     2 
0.88546       0.86708 
OT0T24T*     tf750~40~3 
0.01838* 

Table  26 

Table   of  means  for  the   significant  effects   in  Table   13, 
geographic  areas   3  and  5  for  milk,   using  Tukey's  hsd  test. 

Per  Cent   Pat 

s-   =   0.00258,      q0> 

D5   =    3.98,      hsd   =  0.01027. 

6 
5 

4 

3 

0.84351 
0.84931 
0.86282 
0.87928 

2 

0.88371 

0.04025* 

O.O3O7I* 

0.01989* 

0.00443 

345 
0.87928        0.86382        0.85300 
3763577*     0.02031*     OTO0~959~ 
0.02628        0.01082* 
0.01546 
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Table  26 

(Continued) 

Table  of  means  for  the  significant  effects  in  Table  13, 
geographic  areas  3  and  5  for  milk,  using  Tukey  s  hsd  test. 

Lactation  Number 

s-  =  0.00283,   q0  05  =  4.16,   hsd  =  O.OII77. 

4 
3 
5 
2 
6 

0.84846 
0.84931 
0.85188 
0.85859 
0.86924 

16        2        5 
0.91050   0.86924   0.85859   0.85188 
0.06204*  3702078*  070~1013   0.00342 
0.06119*  0.01993*   0.00928   0.00257 
0.05862*  O.OI736*  0.00671 
0.05191*   0.01065 
0.04126 

3 
0.84931 
0.00085 

Season 

s-  =  0.00251,    qQ  Q5  =  5.74,   hsd  =  0.00864. 

3 

2 
4 

14         2 
0.88474   0.87699   0.85542 
0.84151  I  0.04525*   0.0'55'4"8"*   OYOI59T* 
0.85542   0.02952*   0.02157* 
0.87699  1  O.OO775 

Table  27 

Table  of  means  for  the  significant  effects  in  Table 
geographic  areas  2  and  4  for  fat,  using  Tukey 's  hsd 

17, 
test. 

Geographic  Area 

s-  =  0.00258,    qQ  Q5  ■  2.89,   hsd  =  0.00746. 

4         2 
0.87484  -  0.86629  =  0.00855* 

Per  Gent  Fat 

s-   =  0.00516,    q0_05  =  5.48,   hsd  =  0.01100. 

2 

3 

4       3 
0.87598   0.87428 
0.86145  1  U701455*  0.01285* 
0.87428  I  0.00170 
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Table  27 
(Continued) 

Table  of  means  for  the  significant  effects  in  Table  17, 
geographic  areas  2  and  4  for  fat,  using  Tukey's  hsd  test. 

Lactation  Number 


0.0044-7, 


3  O.85251 
6  0.8554-1 
5  0.85570 

4  0.85606 
2  0.86578 


0.93789 

0.08248* 
0.08219* 
0.08183* 
0.07211* 


0.05 

2 

0.86578 
"OTOT327 
O.OIO37 
0.01008 
0.00972 


=  4.30,   hsd  =  0.01922. 


4 
0.85606 

0.00065 
0.00036 


5 

O.8557O 
O.OO319 
0.00029 


Season 
s-   =  0.00365,   q0-05  =  3-84, 


2  0.86074 
1  0.86810 
4  0.87493 


3 
0.87847 
"0Tor773* 
O.OIO37 
0.00354 


4 
0.87493 

'o.'oi4rg* 

0.0068 3 


hsd   =   0.01402. 
1 


0.86810 
CT0W3S 


Table  28 

Table  of  means  for  the  significant  effects  in  Table  18, 
geographic  areas  3  and  5  for  fat,  using  Tukey's  hsd  test. 

Per  Cent  Fat 


0.00326, 


0.85432 
0.85694 
0.86233 
0.86323 


40.05 


3-< 


5 

0.87182 

O.OI75O* 

0.01488* 

0.00949 

0.00859 


4 
0.86325 
0.0089T 
0.00629 
0.00090 


hsd  =  0.01297. 


6 
0.86233 
0.00801 
0.00539 


0.85694 
0.00262 
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Table_28 
(Continued) 

Table  of  means  for  the  significant  effects  in  Table  18, 
geographic  areas  3  and  5  f°r  fat,  using  Tukey's  hsd  test. 


s-  =  0.00357, 


Lactation  Number 
q0.05  =  4-16>   hsd 


0.01485 


1        2 

0.92421  0.86179 

4  0.84130   "0758251*  0.02049" 
3  0.84297   0.08124*  0.01882*   O.OI338* 

5  0.84374   0.08047*  0.01805*   0.01261 

6  0.85635   0.06786*  0.00544 
2  0.86179   0.06242* 


6         5 
0.85635   0.84374 
5751 515*  0755244" 
0.00077 


s-   =  0.00292, 


2  0.85658 
1  0.85811 
4  0.86531 


3 
0.86691 
O.OIO33 
0.00880 
0.00160 


Season 

90.05  =  3'74' 
4 

0.86531 
0. OO87 3 
0.00720 


hsd  =  0.01092. 


1 
0.85811 
0.00153 


0.84297 
0.00167 


Beginning  with  Per  Cent  Pat  for  milk  in  areas  2  and  4, 
[Table  25] ,  it  can  be  seen  that  levels  2  and  3  differ  signifi- 
cantly from  level  4,  but  not  from  each  other.  This  would  suggest 
that  the  lactation  estimation  curves  for  milk  in  areas  2  and  4 
for  low  fat  producing  Holsteins  and  Holsteins  and  Milking  Short- 
horns are  essentially  the  same,  but  that  they  differ  from  the 
average  curve  for  the  other  breeds.  This  relationship  is  given, 
diagramatically,  in  Pigure  2  below. 

Figure  2 

A  diagram  of  the  mean  comparisons  for  milk  for  Per  Cent  Pat  in 
areas  2  and  4. 


23. 
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This  would  suggest  using  only  two  levels  of  this  factor,  i.e., 

one  for  high  producing  cows  [4.1%  and  over]  and  one  for  low 

producing  cows  [4.0%  and  under]. 

The  analysis  of  means  for  Lactation  Number  in  areas  2  and 

4  for  milk  substantiate  the  results  found  by  Madden  et^  al.  [1956] 

and  other  workers  cited  in  Section  1.2.4,  namely,  that  factors 

for  milk  can  he  separated  into  two  groups,  i.e.,  first  lactations 

and  second  lactations  or  greater.   This  point  is  illustrated  in 

the  following  diagram. 

Figure  3 

A  diagram  of  the  mean  comparisons  for  milk  for  Lactation  Number 
in  areas  2  and  4. 
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The  analysis  of  means  for  Season  indicates  that  curves  for 
cows  freshening  in  the  spring  and  summer  can  be  combined,  giving 
three  groups,  one  for  winter,  one  for  spring  and  summer,  and 
one  for  fall.   Diagramatically,  this  relationship  is  given  in 
Figure  4. 

Figure  4 

A  diagram  of  the  mean  comparisons  for  milk  for  Season  in  areas 

2  and  4. 


Cows  freshening  in  winter  have  a  steeper  lactation  estimation 
curve  than  those  freshening  during  other  seasons  of  the  year, 
and  cows  freshening  in  cold  seasons  have  a  steeper  lactation 
estimation  curve  than  do  those  freshening  in  hot  months.   [See 
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figure  14  and  the  discussion  following] . 

In  the  analysis  of  means  for  areas  3  and  5  for  milk  [Table 

26],  for  Per  Cent  Pat,  levels  2  and  3  and  levels  5  and  6  could 

be  combined,  but  the  other  levels  differ  significantly  from 

each  other.  This  can  be  seen  in  the  following  diagram. 

Figure  3 

A  diagram  of  the  mean  comparisons  for  milk  for  Per  Cent  Pat 
in  areas  3  and  $• 


In  areas  3  and  5,  the  lactation  estimation  curves  for  milk 

could  be  considered  different  for  breeds  corresponding  to  4.0% 

and  under,  4.1%  to  4.6%,  and  4.7%  and  over,  respectively.   In 

areas  2  and  4  cows  producing  4.7%  to  5-2%  combined  with  the 

higher  fat  producing  cows  and  in  areas  3  and  5  they  are  not. 

In  the  former  case,  results  similar  to  those  for  areas  3  and  5 

would  probably  have  been  obtained  also,  had  the  4.7%  to  5-2% 

group  been  kept  separate. 

Por  Lactation  Number,  a  more  complicated  relationship 

is  found  [see  Pigure  13].  Here,  first  lactations  differ  from 

all  other  lactations,  but  curves  for  sixth  lactations  or  over 

also  differ  significantly  from  those  for  all  other  lactations 

[with  the  exception  of  the  second  lactations]. 

Pigure  6 

A  diagram  of  the  mean  comparisons  for  milk  for  Lactation  Number 
in  areas  3  and  5- 
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In  this  case,  the  curves  could  probably  be  grouped  into  first 

lactations,  adult  cow  lactations,  and  old  cow  lactations.   Or 

they  could  be  grouped  so  that  the  second  and  sixth  lactations 

are  estimated  by  the  same  factor. 

A  different  relationship  is  found  for  different  seasons 

in  areas  3  and  5  for  milk.   The  analysis  of  means  indicates 

that  curves  for  cows  freshening  in  the  fall  and  winter  could 

be  combined,  but  that  they  differ  for  cows  freshening  in  the 

spring  and  summer.  See  the  diagram  below. 

Figure  7 

A  diagram  of  the  mean  comparisons  for  milk  for  Season  in  areas 
3  and  5- 


It  can  be  seen  that  cows  freshening  in  the  cold  seasons  have  a 
steeper  lactation  estimation  curve  than  those  freshening  in  the 
hot  seasons,  with  the  steepest  curve  being  for  cows  freshening 
in  the  winter  months.   [See  Figure  14]. 

Cows  from  New  Mexico  and  Arizona  have  steeper  curves  than 
do  those  from  the  northern  states  [Table  27] . 

For  Per  Cent  Fat,  for  fat  an  inverse  relationship  to  that 
noted  for  milk  in  areas  2  and  4  is  found.  The  Holsteins  and 
the  red  breeds  should  be  grouped  together  and  the  low  fat 
producers  kept  separate.  Again  a  diagram  of  the  mean  comparisons 
is  given. 

Figure  8 

A  diagram  of  the  mean  comparisons  for  fat  for  Per  Cent  Fat  in 
areas  2  and  4. 
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The  analysis  of  means  for  Lactation  Number  gives  the  same 
result,  in  terms  of  significance,  for  fat  as  for  milk,  although, 
the  rankings  of  the  means  are  different.  Here  again,  one  need 
only  have  a  factor  for  the  first  lactations  and  one  for  all 
succeeding  lactations. 

A  comparison  of  the  means  for  Season  in  areas  2  and  4-  for 
fat,  indicates  a  somewhat  more  complicated  situation  than  that 
for  milk;  the  following  diagram  may  serve  to  illustrate  the 
situation. 

Figure  9 

A  diagram  of  the  mean  comparisons  for  fat  for  Season  in  areas 

2  and  4. 


This  diagram  of  the  analysis  of  means  indicates  that  lactation 

estimation  curves  for  summer,  fall,  and  winter  do  not  differ 

from  each  other,  and  that  the  curves  for  winter  and  spring  do 

not  differ  from  each  other.   It  would  probably  be  wisest  to 

combine  summer  and  fall,  and  to  combine  winter  and  spring.   It 

should  be  noted  "chat  the  steepest  curve  is  for  cows  freshening 

in  the  summer,  and  the  flattest  curve  is  for  cows  freshening 

in  the  spring  [in  terms  of  estimating  fat  production]. 

Table  28,  which  gives  the  comparisons  of  means  for  fat 

for  areas  J  and  5  yields  the  following  diagram  for  Per  Cent  Fat. 

Figure  10 

A  diagram  of  the  mean  comparisons  for  fat  for  Per  Cent  Fat  in 
areas  3  and  5- 


4  6 


56 


This  diagram  of  the  analysis  of  means  indicates  that  while  the 

average  curve  for  4.7%  to  5-2%  fat  producers  differs  from  those 

for  3-4-%  and  under  and  J. 5%  to  4.0%  fat  producers,  the  curve 

for  4.7%  to  5.2%  does  not  differ  from  those  of  5.3%  and  over 

and  4.1%  to  4.6%  fat  producers,  nor  do  these  latter  two  differ 

from  those  of  3.4%  and  under  and  3.5%  to  4.0%  fat  producers. 

This  abrupt  departure  from  the  linear  trend  which  existed  in 

the  milk  means  as  the  Per  Cent  Fat  increased,  strengthens  the 

idea  that  breed  differences  rather  than  average  per  cent  fat 

affect  the  shape  of  the  lactation  estimation  curve.   It  seems 

more  likely  that  if  the  effect  were  due  to  level  of  per  cent 

fat,  the  trend  would  be  more  linear  [or  slightly  curvelinear  in 

nature],  and  not  have  the  tremendous  drop  in  the  sixth  level, 

which  was  encountered. 

The  relationship  of  the  means  with  respect  to  Lactation 

Number  seems  somewhat  different  in  the  fat  coefficients  than  in 

the  milk  coefficients  [Figure  11] . 

Figure  11 

A  diagram  of  the  mean  comparisons  for  fat  for  Lactation  Number 
in  areas  3  and  5« 

1  2 6     5  34 


This  diagram  brings  out  the  fact  that  if  the  means  are  combined 
as  for  milk,  the  groupings  would  be  the  same,  with  one  factor 
being  used  to  estimate  the  first  lactation  production,  one 
for  the  second  and  sixth  or  greater  lactations,  and  one  for 
the  intermediate  lactations.  The  only  differences  between 
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these  two  tables  [Table  26  and  Table  28]  are  that  the  rankings 
of  the  means  for  the  second  and  sixth  lactations  are  reversed, 
and  the  difference  between  the  second  and  fourth  lactation 
means  is  not  significant  for  milk  as  it  is  for  fat.   Thus,  the 
same  pattern  of  estimation  for  milk  and  fat  could  be  used  in 
structuring  the  factors. 

The  analysis  of  means  for  Season  did  not  yield  any 
significant  differences,  however,  here  as  in  areas  2  and  4, 
curves  for  summer  and  fall  seem  steeper  than  for  winter  and 
spring. 

Graphs  of  the  various  levels  will  now  be  given  as  a 
further  aid  in  the  study  of  the  means  of  significant  effects. 
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The  graphs  in  Figures  12  through  16  show  the  effects  of 
moving  from  one  level  of  each  factor  to  another  in  terms  of  the 
average  or  marginal  factors.  On  each  of  these  graphs,  a  high 
coefficient  indicates  a  steep  lactation  estimation  curve  which 
means  that  the  lactation  curve  for  the  cow  is  quite  flat.  This 
one  coefficient  characterizes  the  entire  lactation  estimation 
curve. 

It  can  he  seen  from  Figure  12  that  the  coefficients  for 
Average  Per  Cent  Fat  for  milk  in  hoth  areas  2  and  4-  and  areas 
3  and  5  show  a  steep  decline  as  one  moves  from  low  fat  producing 
cows  to  high  producing  cows.  This  would  indicate  that  low  fat 
producing  breeds  produce  proportionately  more  milk  later  in 
the  lactation  than  they  do  in  the  early  part;  the  opposite  is 
true  for  the  high  fat  producing  breeds.   As  would  he  expected, 
coefficients  for  these  two  sets  of  areas  for  fat  indicate  that 
as  Average  Per  Cent  Fat  increases,  the  coefficients  also  increase. 
This  indicates  that  lactation  estimation  curves  are  steeper  in 
high  fat  producing  breeds  than  in  low  fat  producing  breeds.  This 
means  that  proportionately  more  fat  is  produced  early  in  the 
lactation  in  low  fat  producing  breeds  than  in  high  fat  producing 
breeds.   Both  Smith  [1959]  and  Johansson  [1961]  have  pointed 
out  that  as  per  cent  fat  increases,  level  of  production  decreases. 
Note  that  the  magnitude  of  this  effect  is  greater  for  milk 
coefficients  than  for  fat  coefficients.  This  indicates  that 
when  per  cent  fat  is  high,  milk  production  is  low,  thus  masking 
the  effects  of  fat  production.  The  sharp  drop  between  the 
fifth  and  sixth  levels  [5.3%  to  5.8%  and  5.9%  and  over]  for 
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fat  in  areas  3  and  5  could  possible  be  indicative  of  a  breed 
difference,  rather  than  a  difference  due  to  change  in  level  of 
per  cent  fat,  as  has  been  mentioned  earlier. 

From  Figure  13,  it  can  be  seen  that  the  effect  of  lactation 
number  is  similar  whether  for  milk  or  fat,  or  whether  in  areas 
2  and  4  or  3  and  5.   From  the  graphs,  there  seems  to  be  a 
relatively  larger  effect  in  areas  3  and  5  than  in  areas  2  and  4. 
This  could  be  due  to  actual  differences  in  the  effect  for  the 
different  areas,  or  it  could  be  due  to  the  fact  that  more 
records  were  used  in  the  calculation  of  means  in  areas  3  and  5 
than  in  areas  2  and  4.  Another  phenomenon  to  be  noted  is  that 
for  milk  coefficients,  the  coefficient  for  the  sixth  and  greater 
lactations  is  consistently  greater  than  that  for  second 
lactations.   For  fat  coefficients,  the  opposite  is  true.   Each 
of  these  curves  seems  to  illustrate  the  fact  that  the  shape  of 
the  lactation  curve  changes  as  a  cow  matures  in  much  the  same 
manner  as  the  total  production  changes.  The  greatest  change 
occurs  between  the  first  and  succeeding  lactations,  the 
steepest  lactation  estimation  curve  being  for  first  lactations. 
The  flattest  lactation  estimation  curves  are  generally  for 
third,  fourth,  and  fifth  lactations.  Except  for  the  case  of 
fat  for  areas  2  and  4,  the  effect  seems  to  be  distinctly 
quadratic  in  nature,  meaning  that  older  cows  have  curves  which 
are  becoming  more  and  more  like  those  of  first  lactation  cows. 
The  effect,  as  expressed  by  the  graph  of  these  coefficients  is 
very  similar  to  the  effect  of  age  on  total  production.  These 
graphs  indicate  that  for  first  lactation  cows,  a  larger 


65 


proportion  of  milk  and  fat  is  produced  later  in  the  lactation 
than  in  the  beginning  portion.  This  proportion  decreases  at  a 
decreasing  rate  until  the  third  and  fourth  lactations,  after 
which  it  again  begins  to  increase. 

Although  the  graphs  shown  in  Figure  14  must  be  interpreted 
in  terms  of  individual  points  for  winter,  spring,  summer,  and 
fall  since  their  consecutiveness  is  not  the  important  fact, 
it  can  be  seen  that  the  graphs  for  milk  are  similar  in  both  sets 
of  areas,  as  are  those  for  fat.   The  coefficients  indicate  that 
curves  for  fail  and  winter  are  similar  and  the  curves  for  spring 
and  summer  are  similar  for  milk,  and  that  curves  for  winter 
and  spring  are  similar  and  those  for  summer  and  fall  are 
similar  for  fat.  This  has  not  been  mentioned  in  earlier  studies. 
In  areas  2  and  4,  the  coefficients  seem  to  be  higher,  in  general, 
for  both  milk  and  fat  than  those  in  areas  3  and  5.   The  co- 
efficient for  milk  is  highest  in  the  winter  and  lowest  in  the 
summer,  and  the  coefficient  for  fat  is  highest  in  the  summer 
and  lowest  in  the  spring. 

In  Figure  15,  it  can  be  seen  that  the  coefficients  for 
milk  follow  a  somewhat  curvelinear  trend  which  increases  as 
one  moves  from  low  to  high  producing  cows.   An  opposite  effect 
is  noted  for  fat,  although  the  effect  is  not  so  pronounced. 
If  a  smooth  curve  is  fitted  through  the  points  on  the  graph, 
it  appears  that  both  effects  are  cubic  in  nature.   The 
coefficients  for  level  of  fat  production  in  Figure  16  exhibit 
the  same  shape  curves  for  milk  and  fat  as  those  in  Figure  15; 
the  same  sort  of  cubic  effect  is  represented.   The  irratic 
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scatter  of  points  is  probably  due  to  random  error. 

5.4.   Results  of  the  analyses  of  interactions  and  discussion. 

Below  are  presented  tables  of  significant  interactions  for  the 

four  four-way  analyses.  Following  these  are  graphs  representing 

each  table. 

Table  29 

Table  of  interactions  for  Geographic  Area  by  Lactation  Number 
for  milk  for  areas  2  and  4.   [See  Table  11.] 

LN 
LEVEL   12       3       4       4       6      AVE 

2  0.92290  0.87363  0.86093  0.87172  0.87319  0.87256  0.87916 

GA  4  0.94573  0.86573  0.86428  0.85412  0.86850  0.87548  0.87897 

AVE  0.93431  0.86968  0.86261  0.86292  0.87085  0.87402  0.87906 

Table  30 

Table  of  interactions  for  Per  Cent  Fat  by  Season  for  milk  for 
areas  2  and  4.   [See  Table  11.] 

S 
LEVEL   12       3       4       AVE 

2  0.89678  0.87199  O.88358  0.89322  0.88639 

■pp  3  0.89939  0.87251  0.87670  0.89589  0.88612 

■"  4  0.90586  0.85673  0.82888  0.86272  0.86468 

AVE  G. 90067  0.86708  0.86305  0.88546  0.87906 

Table  31 

Table  of  interactions  for  Geographic  Area  by  Per  Cent  Fat  for 
milk  for  areas  3  and  5.   [See  Table  12.] 

PF 
LEVEL   2       3       4       5       6       AVE 

3  0.88824  O.88375  0.86734  0.85739  0.83742  0. 86683 

GA  5  0.87918  0.87480  O.86O3O  0.84860  0.84960  0.86250 

AVE  0.88371  0.87928  0.86383  0.85300  0.84351  0.86466 

67 


Table  52 

Table  of  interactions  for  Per  Cent  Fat  by  Lactation  Number  for 
milk  for  areas  3  and  5-   [See  Table  12] . 

LN 
LEVEL   12       3       4       5       6      AVE 

2  0.94251  0.87417  0.86616  0.86626  0.86989  0. 88528  O.8837I 

3  0.93164  0.87512  0.86145  0.86271  0.87011  0.87666  0.87928 
™,  4  0.90926  0.86456  0.85285  0.84800  0.84406  0.86420  0.86382 
r  5  0.89621  0.84142  0.83506  0.85029  0.84436  0.85064  0.85300 

6  0.87288  0.83965  0.83104  0.81506  0.83098  0.87143  0.84351 
AVE  0.91050  0.85859  0.84931  0.84846  0.85188  0.86924  0.86466 


Table  55 

Table  of  interactions  for  Per  Cent  Pat  by  Season  for  milk  for 
areas  3  and  5.   [See  Table  12]. 

S 
LEVEL   12       3       4       AVE 

2  0.88887  0.87524  0.87476  0.89598  O.8837I 

3  0.89578  0.87418  0.86076  0.88641  0.87928 
pp  4  0.88805  0.84752  0.83640  0.88332  0.87928 
r*    5  0.88065  0.84477  0.82669  0.85988  O.853OO 

6  O.87O35  0.83540  0.80893  0.85935  0.84351 
AVE  0.88474  0.85542  0.84151  0.87699  0.86466 


Table  54 

Table  of  interactions  for  Per  Cent  Fat  by  Season  for  fat  for 
areas  2  and  4.   [See  Table  16]. 

S 
LEVEL   12       3       4       AVE 

2  0.85370  0.83871  0.88133  0.87197  0.86143 

3  0.86562  0.86465  O.8857O  0.88115  0.37428 

4  0.88500  0.87889  O.86837  0.87166  0.87598 
AVE  0.86810  O.87074  0.87847  0.87493  O.87056 


Table  35 
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Table  of  interactions  for  Geographic  Area  by  Per  Cent  Fat  for 
fat  for  areas  3  and  5.   [See  Table  17]. 

PF 
LEVEL   2       3       4       5       6       AVE 

3  0.85713  0.86424  0.87297  0.87319  0.85098  0.86370 

GA  5  0.85152  0.84965  0.85350  0.87044  0.87367  0.85976 

AVE  0.85432  0.85694  0.86323  0.87182  0.86233  0.86173 
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The  above  graphs  and  tables  point  out  the  nature  of  the 
interactions.  A  study  of  Figure  17  indicates  that  certain 
differences  do  exist  between  the  behavior  of  Lactation  Number 
in  Area  2  and  that  in  Area  4.  The  difference  between  first 
lactations  in  the  two  areas  indicates  that  Hew  Mexico-Arizona 
cows  have  a  steeper  lactation  estimation  curve  than  those  in 
Area  2,  but  in  terms  of  interaction,  it  is  important  to  note 
that  the  ranks  of  the  coefficients  for  the  two  areas  reverse 
at  every  level  of  Lactation  number  except  one.  The  steepness 
of  the  lactation  estimation  curve  is  a  minimum  at  the  third 
lactation  for  Area  2  and  at  the  fourth  lactation  for  Area  4. 
Also,  the  over-all  effect  of  Lactation  Number  is  greater  in 
Area  4  than  in  Area  2. 

Prom  Figure  18,  one  can  see  a  very  definite  difference 
between  the  effects  of  level  of  Average  Per  Cent  Fat,  depending 
on  the  season  in  which  the  cow  freshens.   For  those  that  freshen 
in  winter  [December-February],  the  higher  the  per  cent  fat,  the 
higher  the  proportion  of  milk  the  cow  produces  later  in  her 
lactation,  although  the  effect  is  not  as  pronounced  as  those 
for  the  other  three  seasons.   This  would  indicate  that  a  cow 
which  produces  a  high  percentage  of  fat  in  her  milk,  freshening 
in  winter,  would  produce  the  largest  proportion  of  her  milk  in 
the  summer  and  fall  months.   A  high  level  of  Per  Cent  Fat  in  a 
cow  freshening  in  the  summer  [June-August]  indicates  that  she 
produces  a  higher  proportion  of  her  milk  during  these  months. 
Coefficients  for  high  fat  cows  are  consistently  lower  than 
for  low  fat  cows.  These  abrupt  differences  seem  especially 
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indicative  of  breed  differences,  rather  than  differences  due 
to  per  cent  fat.   In  Figure  19,  this  difference  in  high  producing 
cows  again  appears  [though  this  time  not  in  terms  of  season]. 
Colorado  [Area  3J  has  a  higher  proportion  of  colored  breeds 
than  does  Utah- Nevada  [Area  5] »  hence,  one  could  say  that  the 
graph  for  GA5  probably  indicates  effects  of  per  cent  fat  change, 
and  the  significant  drop  between  PF5  [4.7-5.2%]  and  PF6  [5.3% 
and  over]  would  indicate  a  breed  difference.  The  other  points 
are  roughly  parallel. 

The  graph  in  Figure  20  seems  to  indicate  that  the  depressing 
effect  on  steepness  of  the  lactation  estimation  curve  caused  by 
an  increase  in  level  of  Per  Gent  Fat  has  much  more  effect  on 
young  cows  than  on  old  ones.  The  abrupt  change  in  trend  dis- 
played by  PF5  and  PF6  could  also  indicate  a  breed  difference, 
as  could  the  lack  of  smoothness  also  noted  in  the  graphs 
between  PF3,  PF4-,  and  PF5. 

The  abrupt  drop  noted  in  Figure  18  does  not  show  up  in 
Figure  21.   This  could  be  because  of  the  fact  that  since  Utah- 
Nevada  is  predominantly  Holstein  what  red-breed  effect  does 
come  from  Colorado  is  masked  by  the  Utah- Nevada  Holsteins. 
It  is  interesting  to  note  the  change  in  SI  between  PF2  and  PF3. 

From  Figure  22,  a  Season  effect  can  be  seen  interacting 
with  Per  Cent  Fat.  As  has  been  earlier  pointed  out,  the  effect 
of  per  Cent  Fat  on  the  milk  lactation  estimation  curve  is 
practically  the  opposite  of  the  effect  on  the  fat  lactation 
estimation  curve.   Instead  of  increase  of  level  of  Per  Cent 
Fat  having  a  depressing  effect  on  the  steepness  of  the 
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lactation  estimation  curve,  it  actually  has  the  opposite  effect. 
Mean  separation  procedures  here  indicated  that  Seasons  1  and  2 
should  be  combined,  as  should  Seasons  3  and  4.  In  this  graph, 
it  can  he  seen  that  these  separations  still  are  indicated.   All 
the  interaction  effect  is  with  respect  to  the  high  fat  cows. 
The  fact  that  Season  interacts  so  strongly  in  this  case  is 
again  indicative  of  the  breed  effect. 

The  effects  of  area  by  Per  Cent  F at  interaction  in  Figure 
23  are  similar  to  those  shown  in  Figure  19.  This  again  is 
probably  indicative  of  the  difference  in  proportions  of  the 
red  breeds  in  Colorado  and  in  Utah-Nevada.   There  could  be  a 
larger  distribution  of  the  red  breeds  in  other  states  in  areas 
2  and  4,  than  in  areas  3  and  5,  and  this  could  probably  account 
for  the  differences  between  graphs  for  interactions  with  respect 
to  Per  Cent  Fat  for  the  two  sets  of  areas,  and  why  Per 'Cent  Fat 
interacted  with  area  for  both  milk  and  fat  only  in  areas  3  and  5. 
3.5  Factors  to  be  used  for  the  lactation  estimation  curve. 
3.5.I.  Preliminaries.  The  discussion  of  the  first  four  sections 
of  this  chapter  indicate  that  a  different  structure  of  factors 
should  be  used  depending  on  the  area  in  which  the  estimates  are 
to  be  used,  and  whether  they  are  for  milk  or  fat.  The  analyses 
have  brought  out  the  fact  that  there  should  definitely  be  a 
different  set  of  factors  for  milk  than  for  fat.  The  discussion 
of  interactions,  all  of  which  were  in  terms  of  Per  Cent  Fat, 
indicate  that  the  fact  that  Breed  is  not  known  tends  to  confuse 
the  effect  of  Area,  and  Area  tends  to  confuse  the  true  inter- 
action of  Breed  with  Season.   It  appears  from  the  analysis  of 
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Geographic  Area  by  Per  Cent  Fat,  that  there  is  probably  an 
effect  due  to  Per  Gent  Pat  and  one  due  to  Breed.  Because  of 
the  lack  of  knowledge  about  Breed  and  its  effect  on  per  cent 
Pat,  and  the  fact  that  Per  Cent  Pat  and  Season  interact  in 
most  cases,  and  that  Per  Cent  Pat  interacts  with  some  of  the 
other  variables,  it  was  decided  to  leave  the  factors  in  their 
original  form  for  use  in  estimation  until  further  work  could 
be  done  with  respect  to  Breed  and  Ber  Cent  Pat.   In  using  cell 
means,  interaction  effects  are  not  confused  in  any  way,  and  the 
estimates  should  not  be  affected  by  the  significant  main  or 
interaction  effects.  The  number  of  factors  is,  of  course, 
somewhat  unweildy.  Factors  will  now  be  given  in  terms  of  each 
of  the  four-four-way  analyses,  and  also  for  the  four  analyses 
in  terms  of  level  of  production. 

3-5.2.   Factors  for  estimating  total  from  part  milk  and  fat 
production  areas  2  and  4.   Tables  25  and  27  and  Figures  2-4,  8-9, 
and  12-14  indicate  that  certain  levels  of  the  variables  could 
possibly  be  combined.  The  differences  in  combination  indicated 
that  a  different  structure  existed  for  the  variables  for  milk 
than  for  fat.  Because  of  the  significant  interactions  involved, 
the  factors  are  not  combined,  but  are  presented  as  cell  means 
for  all  effects. 
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Table  36 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  2  and  4. 

Area  2 

[Worth  Dakota,  Iowa,  Nebraska,  Montana,  Idaho,  and  Wyoming;] 

Average  Per  Cent  Fat- 3.4-%  and  under 

First  Lactations  Milk            Fat 

December-February  0.95141  0.90529 

March-May  0.93451  0.90576 

1  June-August  0.93231  0.9*179 

September-November  0.9*271  0.92937 


Second  Lactations 

December-February  0.87682 

March-May  0.88493 

June-August  O.8743I 

September- November  0.86054 

Third  Lactations 

December-February  0.85978 

March-May  0.83907 

June- August  O.873OO 

September-November  0.87396 

Fourth  Lactations 

December-February  0.9*249 

March-May  0.87040 

June-August  0.86692 

September- November  O.90277 

Fifth  Lactations 

December-February  0.85926 

March-May  0 . 88440 

June-August  0.92110 

September- November  O.878I7 

Sixth  Lactations  and  over 

December-February  0.89862 

March-May  0.82306 

June-August  0.87520 

September-November  0.90789 


0.83088 
0.85752 
O.903OO 
O.83723 


0.81853 
0.87518 
0.87518 
0.84391 


0.88141 
0.81032 
0.89288 
0.88121 


0.81845 
0.85394 
O.87I8O 
0.85891 


0.85227 
0.78282 
0.86582 
0.86650 


First  Lactations 

December-February 
March-May 
June- August 
September-November 


Average  Per  Cent  Fat-3.5%  to  4.0% 


0.94585 
0.89650 
0.90769 
0.94451 


0. 90855 
0.88448 
0.93700 
0.93761 
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Table  36 

[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  2  and  4. 

Area  2 
[Worth  Dakota,  Iowa,  Nebraska,  Montana,  Idaho,  and  Wyoming] 

Average  Per  Cent  Fat- 5. 5%  to  4.0% 

Second  Lactations  Milk            Fat 

December-February  0.90729  0.88052 

March-May  0.88557  0.86101 

June-August  0.87604  0.90006 

September-November  0.88956  0.87354 

Third  Lactations 


December-February 

0.86528 

0.82374 

March-May 

0.85938 

0.809 34 

June-August 

0.84846 

0.86863 

September- November 

0.89517 

0.88255 

Fourth  Lactations 

December- February 

0.87483 

0.82242 

March-May 

0.86661 

0.82709 

June-August 

0.83494 

0.82746 

September-November 

0.88084 

0.85975 

Fifth  Lactations 

December-February 

0.90930 

0.84557 

March-May 

0.83704 

0.81281 

June- August 

0.86874 

0.87902 

September-November 

0.96311 

O.83393 

Sixth  Lactations  and  over 

• 

December-February 

0.90607 

0.82480 

March-May 

0.86636 

0.86008 

June-August 

0.90036 

0.90770 

September- November 

0.90241 

0.87596 

Average  Per  Cent 

Fat-4.1%  and  over 

First  Lactations 

December-February 

0.94890 

0.94258 

March-May 

0.89103 

0.93601 

June-August 

0.88952 

0.93521 

Sept ember- November 

0.88982 

0.91253 

Second  Lactations 

December-February 

0.88650 

0.84565 

March-May 

0.85278 

0.83405 

June-August 

0.82561 

0.88423 

Sept ember- November 

0.86586 

0.86406 
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Table  36 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  2  and  4. 

Area  2 
[North  Dakota,  Iowa,  Nebraska,  Montana,  Idaho,  and  Wyoming] 

Average  Per  Cent  Fat-4.1%  and  over 

Third  Lactations                 Milk  Fat 

December-February         0.87980  0.8 3808 

March-May                0.86639  0.87266 

June-August               0.84085  0.90124 

September- November        0.82609  0.83080 

Fourth  Lactations 


December-February 

March-May 

June-August 

Sept ember- November 

Fifth  Lactations 

December- February 
March-May 
June-August 
September- November 

Sixth  Lactations  and  over 
December-February 
March-May 
June-August 
Sept ember- November 


0.92603 
0.82359 
0.83214 
0.8 3908 


0.88592 
0.87697 
0.80288 
0.89142 


0.96896 
O.8O575 
0.7384-2 
0.87766 


0.87225 
0.88061 
0.864-96 
0.84-230 


0.85863 
0.92573 
0.02606 
0.90557 


0. 90188 
0.82959 
0. 74-906 
0.87135 


Area  4 
[New  Mexico  and"  Arizona] 


First  Lactations 

December-February 
March-May 
June-August 
September- November 

Second  Imitations 

December-February 
March-May 
June- August 
September- November 

Third  Lactations 

December-February 
March-May 
June-August 
September- November 


Average  Per  Cent  Fat-5.4%  and  under 


0.97277 
0.93617 
0.95232 
0.97051 


O.87471 
0.86168 
0.86758 
0.89332 


0.88339 
0.85138 
0.84021 
0.86386 


0.94521 
0.92246 
0.97166 
0.96907 


0.83975 
0.83855 
0.87342 
0.88776 


0.83206 
O.83793 
0.83022 
0.83616 
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Table  56 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  2  and  4. 

Area  4 


[New  Mexico  and  Arizona] 

Average 

Per  Cent  Fat- 3. 

.4%  and  under 

Fourth  Lactations 

Milk 
0.87092 
0.84171 
0.81703 
0.84191 

Fat 

December-February 
March-May 
June-August 
September- November 

0.85025 
0.79888 
0.79820 
0.84985 

Fifth  lactations 

December-February 
March-May 
June- August 
September- November 

0.87974 
0.89097 
0.36508 
0.87646 

0.85815 
0.82997 
0.85854 
0.85720 

Sixth  Lactations  and  over 
December-February 
March-May 
June-August 
September- November 

0.89140 
0.84556 
0.91791 
0.90110 

0.85214 
0.81064 
0.91371 
0.86652 

Averap 

;e  Per  Cent  Fat- 

•  3.5%  to  4.0% 

First  Lactations 

December- February 
March-May 
June-August 
September- November 

0.97050 
0.94917 
0.95099 
0.96816 

0.96472 
0.94915 
0.96786 
0.88562 

Second  Lactations 

December-February 
March-May 
June -August 
September-November 

0.88920 
0.86380 
0.85449 
0.88341 

0.88598 
0.87910 
0.86276 
0.88562 

Third  Lactations 

December-February 
March-May 
June-August 
September- November 

0.86550 
0.84865 
0.87421 
0.88856 

0.85454 
0.84526 
0.87608 
0.87372 

Fourth  Lactations 

December-February 
March-May 
June-August 
September- November 

0.91909 
0.86901 
0.86802 
0.85949 

0.89125 
0.89048 
0.87390 
0.84202 
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Table 

36 

[continued] 

Factors  for  use  in  estimating 

total  from  part  milk  and  fat 

production  in  areas  2  and  4. 

Area 

4 

[New  Mexico  and  Arizona] 

Average 

Per  Cent  Fat- 3. 5% 

to  4.0% 

Fifth  Lactations 

Milk 
0.87408 

Fat 

December-February 

0. 84922' 

March-May 

0.85918 

0.86283 

June-August 

0.86391 

0.85504 

September- November 

0.88228 

0.84904 

Sixth  Lactations  and  over 

December-February 

0.86984 

0.83272 

March-May 

0.86839 

0.89593 

June-August 

0.87251 

0.87294 

Sept ember- November 

0.98333 

0.88265 

Average 

Per  Cent  Fat-4.1% 

and  over 

First  Lactations 

December-February 

0.93383 

0.94384 

March-May 

0.90562 

0.92525 

June-August 

0.89666 

0.93210 

September- November 

0.94210 

0.96234 

Second  Lactations 

December-February 

0.88480 

0.89515 

March-May 

0.84691 

0.85428 

June- August 

0.83652 

0.86191 

Sept ember- November 

O.83235 

0.84498 

Third  Lactations 

December-February 

0.89360 

0.90485 

March-May 

0.85200 

0.85920 

June-August 

0.85356 

0.87055 

September-November 

0.85870 

0.85799 

Fourth  Lactations 

December-February 

0.86531 

0.86761 

March-May 

0.82002 

0.87342 

June-August 

0.81368 

0.90299 

September-November 

0.86325 

0.86401 

Fifth  Lactations 

December-February 

0.89445 

0.87111 

March-May 

0.90703 

0.92643 

June-August 

0.76627 

0.82402 

September-November 

0.86257 

0.86510 
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Table  56 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  2  and  4. 

Area  4 

[Mew  Mexico  and  Arizona] 

Average  Per  Cent  Fat-4.1%  and  over 

Sixth  Lactations  and  over  Milk            Fat 

December-February  0.90226  0. 87833 

March-May  0.83264  0.82946 

June-August  0.85242  0.86811 

September-November  0. 85836  0.83891 


3.5.3.   Factors  for  estimating  total  from  part  milk  and  fat 
production  in  areas  3  and  5*  Tables  26  and  28  and  Figures 
5-7,  10-11,  and  12-14  indicate  that  certain  levels  of  the 
variables  could  be  combined.  The  differences  in  combinations 
again  indicate  basic  differences  in  structure  between  milk  and 
fat.  The  differences  between  areas  3  and  5  is  probably  due 
largely  to  the  difference  in  relative  proportion  of  breeds 
between  the  two  areas  [see  Section  3.4-3  and  to  interaction  of 
Breed  and  Per  Cent  Fat  with  other  variables,  particularly 
Season.   Because  of  the  significant  interactions  involved,  the 
factors  are  not  combined,  but  are  presented  as  cell  means  for 
all  effects. 
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Table  57 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  3  and  5- 

Area  3 

[Colorado] 

Average  Per  Cent  Fat-  3-4-%  and  under 

Hi  Ik  • Fat 


First  Lactations 

December-February  0.94801 

March-Hay  0 . 9  3690 

June-August  0.94314 

September-November  0.97687 

Second  Lactations 

December-February  0.88845 

March-May  0. 87234 

June-August  0.86611 

September- November  0.89103 

Third  Lactations 

December-February  0.87545 

March-Hay  0.87797 

June-August  0.86045 

September- November  0.88624 

Fourth  Lactations 

December-February  0.87504 

Harch-Hay  0.84229 

June-August  0.86304 

September-November  O.87OIO 

Fifth  Lactations 

Decembef^February  0.87229 

March-May  0. 85258 

June-August  0.87151 

September- November  0.87546 

Sixth  Lactations  and  over 

December-February  0.87002 

March-May  0.89227 

June-August  0 . 90886 

September- November  0.90143 


.90233 
0.91171 
0.93929 
0.94981 


0.85541 
0.84906 
0.87218 
0.85482 


0.83142 
0.84072 
0.84365 
0.86172 


0.81749 
0.81719 
0.85181 
0.82145 


0.83309 
O.83257 
0.84175 
0.82620 


0.81786 
0.84888 
0.89634 
O.85432 


First  Lactations 

December-Febru ary 
March-Hay 
June-August 
September- November 


Average  Per  Cent  ffat-5.5%  to  4.0% 


0.94415 
0.92286 
0.93070 
0.95662 


0.91706 
0.91704 
0.94650 
0.93505 
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Tattle  37 
Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  3  and  5- 

Area  3 
[Colorado] 

Average  Per  Cent  Fat- 3.%  to  4.0% 


Second  Lactations 

Decemtter-FeEruary 
March-May 
June-August 
September-November 

Third  Lactations 

Decemtter-Fettruary 
March-May 
June- August 
Septemtter-November 

Fourth  Lactations 

December-February 
March-May 
June-August 
September-November 

Fifth  Lactations 

Decemtter-Fettruary 
March-May 
June-August 
September-November 

Sixth  Lactations  and  over 
Decemtter-Fettruary 
March-May 
June-August 
Septemtter-Novemtter 

First  Lactations 

December-Fettruary 
March-May 
June-August 
September- November 

Second  Lactations 

December-February 
March-May 
June-August 
September- November 


Milk 
0.89992 
0.86982 
0.8574-7 
0.894-09 


0.87155 
0.86726 
0.84184 
0.8704-7 


0.90326 
0.87015 
0.84392 
0.85579 


0.89150 
0.88049 
0.84852 
0.87878 


0.88858 
0.86636 
0.86814 
0.88790 


.Fat 
0.86836 
0.85799 
0.86237 
O.87832 


0.82815 
0.85267 
0.83772 
0.83296 


0.85692 
0.84516 
0.83993 
0.83736 


O.85435 
0.85952 
0.84518 
0.84243 


0.85327 
0.87237 
0.85770 
0.84342 


Average  Per  Cent  Fat-4.1%  to  4.6% 


0.94226 
0.90472 
0.87944 
0.93037 


0.89928 
0.86949 
0.84575 
0.88078 


0. 93816 
0.92676 
0.91615 
O.94443 


0.89063 
0.88939 
0.88978 
0.88362 


Table  57 
Factors  for  use  in  estimating  total  from  part  milk  ana  fat 
production  in  areas  3  a*icl  5« 

Area  5 
[Colorado] 


Average  Per  Cent  Fat-4.1%  to  4.6% 
Third  Lactations  Milk  ggt 

December-February         0.88063         0.- 85789 

March-May  0.82383  0.82825 

June-August  0.82909         0.87589 

September-November         0.87916         0.86218 

Fourth  Lactations 

December-Fe  bruary 
March-May 
June -August 
September-November 

Fifth  Lactations' 

December-February 
March-May 
June-August 
September- November 

Sixth  Lactations  and  over 
December-February 
March-May 
June-August 
September- November 

First  Lactations 

December-February 

March-May 

June-August 

Sept ember- November 

Second  Lactations 

December-February 
March-May 
June-August 
September- November 

Third  Lactations 

December-February 
March-May 
June-August 
September- November 


0.88627 
0.84186 
0.79914 
0.85912 

0.86109 
0.86211 
0.79057 
0.84817 

0.87543 
0.81644 
0.80956 
0.88040 

0.84885 
0.85177 
0.86513 
0.87748 

0.88557 
0.84590 
0.85201 
0.89990 

0.85081 
0.85953 
0.86224 
0.894-54- 

Average  Per  Cent 

Fat-4.7%  to  5.2% 

0.94903 
0.87499 
0.88451 
0.91333 

0.93624 
0.89358 
0.95502 
0.94-976 

0.87136 
0.84463 
0.83785 
0.86365 

0. 861 54- 
0.87035 
0.87729 
0.89091 

0.87741 
0.85978 
0.77404 
0.87550 

0.85538 
0.83574- 
0.82743 
0.88111 
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Table  37 
[continued] 

Factors  for  use  in  estimating 
production  in  areas  3  and-  5- 

total  from  part 

milk  and  fat 

Area  3 
[Colorado] 

Average 

!  Per  Cent  Fat-4, 

,7%   to  5-2% 

Fourth  Lactations 

December-February 
March-May 
June-August 
September- November 

Milk 
0.88909 
0.86048 
0.84025 
0.84427 

•Fat 

0".'88'08'0" 
0.86232 
0.91914 
0.85796 

Fifth  Lactations 

December-February 
March-May 
June-August 
September-November 

0.86703 
0.82014 
0.79757 
0.82974 

O.83433 
0.86371 
0.81 305 
0.82406 

Sixth  Lactations  and  over 
December-February 
March-May 
June-August 
September- November 

0.87770 
0.83520 
0.88240 
0.82938 

0.87714 
0.83091 
0.91964 
0.83946 

First  Lactations 

December-February 
March-May 
June-August 
September- November 

Avers 

ige  Per  Cent  Fat- 

0.90768 
0.85092 
0.80793 
0.87951 

■5-3%  and  over 

0.90295 
0.85176 
0.91931 
0.90337 

Second  Lactations 

December-February 
March-May 
June-August 
September-November 

0.87896 
0.82964 
0.7966^ 
0.85215 

0.86646 
O.83O87 
0.85690 
0.85429 

Third  Lactations 

December-February 
March-May 
June-August 
September- November 

0.85910 
0.78826 
0.80436 
0.83255 

0.84919 
0.80991 
0.86272 
0.84949 

. 

Fourth  Lactations 

December-February 
March-May 
June-August 
September-November 

0.809 30 
0.85080 
0.74595 
0.83255 

0.81816 
0.89435 
0.74782 
0.83553 
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Table  57 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  3  and  5- 

Area  5 
[Colorado] 


Average  Per  Cent  Fat-5-3%  and  over 


Fifth.  Lactations 

December-February 
March-May 
June-August 
September-November 

Sixth  Lactations  and  over 
December-February 
March-May 
June- August 
Sept ember- November 


Milk 
0.81737 
0.84924 
0.81101 
0.83645 


0.87260 
0.84191 
0.85028 
0.90363 


Fat 
0.7y804 
0.89982 
0.84199 
0.80185 


0.86680 
0.83981 
0.84056 
0.88162 


Area  5 
[Utah  and  .Nevada] 


First  Lactations 

December-February 

March-May 

June- August 

Sept ember- November 

Second  Lactations 

December-February 
March-May 
June-August 
September- November 

Third  Lactations 

December-February 
March-May 
June-August 
September- November 

Fourth  Lactations 

December-February 
March-May 
June-August 
September-November 


Average  Per  Cent  Fat- 3. 4%  and  under 


0.94848 
0.92042 
0.91874 
0.94746 


0.87976 
0.85985 
0.85817 
0.88766 


0.87075 
0.84668 
0.83825 
0.87346 


0.87173 
0.86776 
0.85893 
0.88121 


0.91579 
0.91416 
0.93013 
0.92504 


0.. 841 33 
0.84084 
0.85834 
0.87006 


0.81994 
0.81092 
0.83417 
0.84741 


0.82786 
0.82983 
O.85154 
0.84049 
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Table  57 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  3  smd  5- 

Area  5 
[Utah  and  Nevada] 


Average  Per  Cent  Fat- 5. 4%  and  under 


Fifth  Lactations 

December-February 
March-May 
June-August 
September- November 

Sixth  Lactations  and  over 
December-February 
March-May 
June-August 
September- November 


First  Lactations 

December-February 
March-May 
June-August 
September-November 

Second  Lactations 

December-February 
March-May 
June-August 
September-November 

Third  Lactations 

December-February 
March-May 
June-August 
September-November 

Fourth  Lactations 

December-February 
March-May 
June-August 
September- November 

Fifth  Lactations 

December-February 
March-May 
June-August 
September-November 


Average 


Milk 

Pat 

0.8860 3 
0.86094 
0.96444 
0.87590 

0.82280 
0.82673 
0.84747 
0.84083 

0.880 36 
0.87290 
0.85553 
0.88491 

0.81812 
0.83547 
0.85617 
0.85101 

i  Per  Cent 

Pat- 5- 5%  to  4.C 

0.93592 
0.91176 
0.91057 
0.94253 

0. 90651 
0.91519 
0.933^-9 
0.93186 

0.87608 
0.85489 
0.84995 
0.88275 

0.84104 
0.85851 
0.85915 
0.85687 

0.87879 
0.85307 
0.85451 
0.87418 

0.82594 
0.82045 
0.84498 
0.85505 

0.88076 
0.85698 
0.84185 
0.84895 

0.82277 
0.82775 
0.85964 
0.80410 

0.88045 
0.85585 
0.85240 
0.87293 

0.80872 
0.82599 
0.84628 
0.82858 

0% 
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Table  37 
[continued] 

Factors  for  use  in  estimating 
production  in  areas  3  ssid.   5- 

total  from  part 

milk  and  fat 

[Utah 

Area  5 

and  Nevada] 

Average  Per  Cent  Fat- 

-5.5%  to  4.0% 

Sixth  Lactations  and  over 
December-February 
March-May 
June-August 
Sept ember- November 

Milk 
0.90037 
0.88065 
0.84943 
0.87188 

'Fat 

0.84609 
0.84674 
0.85522 
0.85290 

First  Lactations 

December-February 
March-May 
June-August 
September- November 

Average  Per  Cent  Fat- 

0.90945 
0.89863 
0.88099 
0.92821 

■4.1%  to  4.6% 

0.89810 
0.92431 
0.90958 
0.92359 

Second  Lactations 

December-February 
March-May 
June-August 
September- November 

0.87754 
0.85577 
0.84643 
0.86259 

0.84885 
0. 79908 
0.87093 
0.85513 

Third  Lactations 

December-February 
March-May 
June-August 
September- November 

0.86541 
0.83577 
0.84821 
O.86O7O 

0.82389 
0. 79908 
0.83405 
0.84585 

Fourth  Lactations 

December-February 
March-May 
June-August 
September- November 

0.86970 
0.82685 
0.83768 
0.86340 

0.83887 
0.85627 
0.85574 
0.84150 

Fifth  Lactations 

December-February 
March-May 
June-August 
September- November 

0.87681 
0.84352 
0.77762 
0.87267 

0.85712 
0.86281 
0.85128 
0.84589 

- 

Sixth  Lactations  and  over 
December-February 
March-May 
June -August 
September- November 

0.88846 
0.82857 
0.83081 
0.88255 

0.84888 
0.84060 
0.85756 
0.84771 
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Table  37 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  3  and  5- 

Area  5 

[Utah  and  Nevada] 

Average  Per  Cent  Fat-4.7%  ,to  3.2% 

Milk  Fat 

TJT&mE  '0  ."92045 

0.89249  0.95497 

0.84644  0.92848 

0.90990  0.94942 


First  Lactations 

December-February 
March-May 
June-August 
September- November 

Second  Lactations 

December-February  0. 83885 

March-May  0.82222 

June-August  0.81029 

September-November  0.84253 

Third  Lactations 

December-February  0. 87432 

March-May  0.81 367 

June- August  0.80241 

September-November  0.82532 

Fourth  Lactations 

December-February  0.86611 

March-May  0 . 80688 

June-August  0.8191 3 

September-November  O.876II 

Fifth  Lactations 

December-February  0 . 86042 

March-May  0.86232 

June-August  0.86235 

September-November  0. 85532 

Sixth  Lactations  and  over 

December-February  0.89754 

March-May  0 . 864  39 

June- August  0.76302 

September-November  0. 85547 

Ave: 
First  Lactations 

December-February  0.89870 

March-May  0.89983 

June-August  0.86601 

September-November  0.85532 


0.84336 
0.85068 
0.87202 
0.861 31 


0.86090 
O.83927 
0.86298 
0.82489 


0.86721 
0.83264 
0.85288 
0.86729 


0.84438 
0.88174 
0.91196 
0.84073 


0.88528 
O.86323 
0.81346 
0.86104 


Average  Per  Cent  Fat- 5. 3%  and  over 


0.91197 
0.94857 
0.94769 
O.923OO 
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Table  57 
[continued] 

Factors  for  use  in  estimating  total  from  part  milk  and  fat 
production  in  areas  3  and  5. 

Area  5 

[Utah  and  Nevada] 

Average  Per  Cent  Fat-5-5%  and  over 

Milk  Fat 


Second  Lactations 

December-February  0.85561 

March-May  0.83225 

June-August  0.81440 

September-November  0.85756 

Third  Lactations 

December-February  0.88712 

March-May  0.82792 

June-August  O.8I5I7 

September-November  0. 83387 

Fourth  Lactations 

December- February  0.87062 

March-May  O.8O378 

June-August  0.78524 

September-November  0.83201 

Fifth  Lactations 

December-February  0. 88325 

March- May  0.80934 

June-August  O.7725O 

September- November  0.86966 

Sixth  Lactations  and  over 

December-February  O.90387 

March-May  0.85095 

June- August  0. 83868 

September- November  0.90954 


0.86477 
0.86051 
0.87190 
0.87734 


0.88287 
0.86089 
0.88322 
0.83992 


0.84704 
0.80724 
0.84819 
0.84819 


0.85886 
0.83069 
0.85461 
0.87328 


0.87715 
0.85654 
0.88072 
0.91289 


3. 5-4.   Factors  for  estimating  total  from  part  production  for 
milk  and  fat  in  terms  of  level  of  milk  and  level  of  fat  production. 
The  factors  here  presented  for  use  in  estimating  total  from  part 
lactation  production  in  terms  of  level  of  milk  and  fat  production, 
although  the  effects  are  significant,  are  probably  of  more 
academic  than  practical  interest.  The  difficulty  in  finding  a 
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base  for  level  of  production  when  the  lactation  is  incomplete 

is  but  one  of  the  problems.  For  this  particular  study,  these 

variables  were  considered  separately.   This  makes  it  difficult 

to  use  them  in  connection  with  the  other  significant  variables. 

It  was  nevertheless  felt  to  be  beneficial  to  include  factors 

for  these  effects  in  a  table  to  make  them  accessible.   Tables 

38  and  39  present  the  factors  for  these  variables. 

Table  38 

Factors  for  use  in  estimating  total  from  part  production  for 
both  milk  and  fat  production  for  level  of  milk  production, 
[in  pounds.] 


Milk  Coefficient 


Fat  Coefficient 


5,999 

and  under 

0.83693 

0.87885 

6,000 

to  6,999 

0.84723 

0.87871 

7,000 

to  7,999 

0.86040 

0.88212 

8,000 

to  8,999 

0.87060 

0.88320 

9,000 

to  9,999 

0.87758 

0.87932 

10,000 

to  10,999 

0.88215 

0.87966 

11,000 

to  11,999 

0.88772 

0.87994 

12,000 

to  12,999 

0.88582 

0.87153 

13,000 

to  13,999 

0.88377 

0.86681 

14,000 

to  14,999 

0.88540 

0.85751 

15,000 

to  15,999 

0.88746 

0.86236 

16,000 

to  16,999 

0.88401 

0.85458 

17,000 

to  17,999 

0.89618 

0.86129 

18,000 

to  18,999 

0.89121 

0.85755 

19,000 

to  19,999 

0.90211 

0.85814 

20,000 

and  over 

0.90526 

0.86579 
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Table  59 

Factors  for  use  in  estimating  total  from  part  production  for 
both  milk  and  fat  production  for  level  of  fat  production 
[in  pounds]. 


309  and  under 

310  to  339 
34-0  to  369 
370  to  399 
400  to  429 
430  to  459 
460  to  489 
490  to  519 
520  to  549 
550  to  579 
580  to  609 
610  to  639 
640  to  669 
670  to  699 
700  and  over 


Milk  Coefficient 

0.86767 
0.87744 
0.87525 
0.88173 
0.88290 
0.88201 
0.88113 
0.88299 
0.88201 
0.88549 
0.89022 
0.88916 
0.88984 
0.90022 
0. 90060 


Fat  Coefficient 

0.88294 
0.881 38 
0.87609 
0.88246 
0.87532 
0.87323 
0.86874 
0.86582 
0.86192 
0.86269 
0.86355 
0.85955 
0.85802 
0.8660 3 
0.85730 


3.6.  Summary  of  the  discussion  of  the  results.  The  analyses 
of  variance  indicated  that  all  variables  studied  had  an  effect 
on  the  lactation  estimation  curve,  either  directly,  or  through 
interaction  with  other  variables,  with  the  exception  of  Days  Dry 
and  Days  Open.  The  significant  interactions  indicated  that  Per 
Cent  Fat  interacted  with  climate  either  by  interacting  with 
Season  of  Freshening  or  by  interacting  with  Geographic  Area, 
although  in  the  latter  case,  the  interaction  could  be  due  to 
breed  differences.   An  analysis  of  the  components  of  variance 
indicated  that  Lactation  Number  accounted  for  the  largest 
proportion  of  the  variation,  in  one  case  accounting  for  more 
than  90%.  Per  Cent  Fat  and  Season  of  Freshening  accounted  for 
relatively  more  variation  for  milk  coefficients  than  for  fat 
coefficients.  An  analysis  of  means  indicated  that  certain 
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variables  could  be  combined,  provided  they  were  involved  in  no 
interaction  effects.   Graphs  of  the  coefficients  indicated  that 
coefficients  for  fat  reacted  similarly  in  both  sets  of  areas, 
as  did  those  for  milk.  While  both  sets  of  coefficients  for  milk 
and  fat  for  lactation  Number  in  areas  2  and  4-  upheld  the  previously 
discussed  results  that  first  lactations  should  be  separated  from 
the  succeeding  ones  and  that  the  succeeding  ones  could  be 
estimated  together,  in  areas  3  and  5,  it  was  found  that  the 
second  and  sixth  or  greater  lactations  should  also  be  estimated 
by  a  separate  coefficient,  although  in  the  case  of  milk,  separate 
estimation  of  the  second  lactation  from  the  third,  fourth,  and 
fifth,  does  not  seem  to  be  quite  so  justified.   It  was  found 
that  a  different  combination  of  seasons  should  be  used  for  milk 
than  for  fat  with  December  through  May  and  June  through  November 
being  used  for  fat  and  September  through  February  and  March 
through  August  being  used  for  fat. 

A  cubic  trend  was  found  to  exist  for  b.oth  Level  of  Milk 
Production  and  Level  of  Fat  Production,  for  both  milk  and  fat 
coefficients.  The  trend  was  in  an  upward  direction  for  milk 
coefficients  for  both  level  of  milk  and  level  of  fat  production, 
and  in  a  downward  direction  for  the  fat  coefficients  in  both 
cases.  The  effect  was  not  as  pronounced  for  Level  of  Fat 
Production  as  for  Level  of  Milk  Production. 

Graphs  of  the  significant  interactions  showed  the  effects 
of  moving  across  all  levels  of  one  of  the  variables  while  holding 
the  other  variable  constant  for  each  of  its  levels.   It  was 
interesting  to  note  the  effects  of  Per  Cent  Fat  with  respect 
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to  interaction.  Evidently  there  was  a  "muddying"  effect  of 
Breed  on  Per  Cent  Fat.  Per  Cent  Pat  interacted  almost 
consistently  with,  season  effects  [the  exception  being  for  fat  in 
areas  3  and  5,  where  there  was  quite  a  large  Per  Cent  Pat  by 
Area  interaction].   Because  of  the  significant  interaction 
effects,  which  influenced  the  results  of  the  mean  comparison 
tests,  one  of  the  variables  was  combined  with  respect  to  level, 
but  the  coefficients  were  left  in  terms  of  cell  means. 

From  all  analyses,  it  was  seen  that  the  structure  of  the 
variables  for  milk  was  distinctly  different  than  that  for  fat. 
The  structural  differences  between  areas  2  and  4  and  areas  3 
and  5  could  possibly  have  been  due  to  the  differences  in  which 
the  levels  were  constructed  for  Per  Cent  Fat  for  the  two  sets 
of  areas,  and  also  to  the  smaller  numbers  of  observations 
involved  in  areas  2  and  4. 

If  one  were  to  use  one  set  of  factors  common  to  the  western 
states,  the  factors  for  Utah  and  Nevada  wou}.d  probably  be  the 
best  to  use,  as  they  are  based  on  the  largest  number  of  obser- 
vations, and  because  the  majority  of  cows  are  of  the  Holstein 
breed. 

As  was  stated  before,  the  factors  for  level  of  production 
are  more  of  an  academic  than  practical  interest  because  of  the 
fact  that  in  their  present  form,  they  cannot  be  used  in 
connection  with  the  other  variables. 
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APPENDIX  I 
On  the  following  pages  will  be  found  examples  of  the 
various  forms  of  the  data  as  it  was  used  as  the  input  and 
obtained  as  output  for  the  principal  programs  and  analyses 
described  in  the  section  on  methods  and  materials.  After 
each  example  will  be  found  a  description  or  discussion  of  the 
data. 
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Figure  24  shows  the  form  of  the  data  after  they  had  been 
loaded  on  tape  and  sorted  into  sequence  hy  year  and  month  within 
cow.   The  record  for  the  cow  which  is  marked  represents  17 
monthly  production  records  beginning  when  the  cow  was  142  days 
into  her  lactation.   The  "H"  is  the  flag  mentioned  in  the  sec- 
tion on  methods  which  indicates  that  either  the  lactation  is 
complete  or  that  the  cow  has  milked  one  month  past  3O5  days. 
For  this  particular  case ,  the  first  part  of  the  lactation  is 
missing,  and  thus,  this  part  of  the  cow's  record  would  not  be 
used.  Also,  the  lactation  was  completed  at  301  days,  hence, 
if  the  early  part  were  present,  the  record  would  be  a  "short 
lactation"  and  still  could  not  be  used.   The  second  lactation 
for  this  cow  also  cannot  be  used,  as  the  latter  portion  of  the 
lactation  is  not  available.  This  is  the  basic  type  of  record 
from  which  the  data  were  taken.   These  records  were  the  input 
to  the  edit  program. 
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Figure  25  is  the  output  to  the  edit  program  and  the  input 
to  the  regression  setup  program  [see  Figure  26].   There  is  a 
great  deal  of  information  in  these  data.  As  can  he  seen  the 
record  hegins  with  the  information  necessary  to  set  up  the 
variable  codes  for  the  8  variables  to  he  studied.   Says  in 
mills  has  been  rolled  back  to  3O5  days  and  a  305-Day  record  has 
been  computed.   The  average  per  cent  fat  is  included  as  a  part 
of  the  305-Day  record.  Following  this  is  given  the  per  cent 
fat  and  pounds  of  milk  produced  for  the  current  period  [the 
milk  is  expressed  to  nearest  10  pounds]  and  the  cumulative 
days  in  milk,  pounds  of  milk,  and  pounds  for  that  period  for 
eight  to  eleven  consecutive  periods.   [The  cumulative  milk 
is  expressed  to  the  nearest  10  pounds,  and  the  cumulative  fat 
to  the  nearest  pound.]  Also  included  is  the  centering  day 
and  the  test  day  for  each  of  the  first  10  periods. 
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Figure  26  is  the  input  to  the  regression  program  mentioned 
on  page  34.   The  data  were  taken  from  the  output  to  the  Data 
Edit  program  and  transformed  into  the  format  shown  above.   In 
doing  so,  the  data  were  coded  according  to  the  8  variables  to 
he  studied.   These  codes  are  found  in  the  area  labeled  "Variable 
Code,"  and  represent  Geographic  Area  [1st  column],  Per  Cent  Fat 
Level  [2nd  column] ,  Lactation  Number  [3rd  column] ,  Season  [4th 
column],  Previous  Days  Dry  [5th  column],  Days  open  [6th  column], 
Level  of  Milk  Production  [7-8th  columns] ,  and  Level  of  Fat 
Production  [9-10th  columns].  The  first  area  is  a  computer 
code  designating  the  type  of  record  and  the  number  of  tape 
positions  in  it.   It  is  for  machine  use  only.   The  third  field 
or  area  is  the  cumulative  days  in  milk  to  that  particular 
period.   The  fourth  field  is  total  305-Day  milk  production. 
The  fifth  is  the  cumulative  milk  production  to  the  given 
period.   The  sixth  is  the  total  305-Day  fat  production;  the 
seventh  is  the  cumulative  fat  production  to  the  given  period. 
The  eighth  is  period  for  which  the  cumulative  record  is  given. 
The  ninth  is  the  total  number  of  periods  in  the  lactation,  or 
305-Day  record.   The  final  field  is  the  cow  number,  a  number 
which  was  assigned  in  consecutive  order  to  the  cows  used  in 
the  study.  The  majority  of  the  305-Day  records  used  had  11 
periods  contained  in  them. 

Note  that  Cow  Number  1  had  3I8  days  in  her  record.  This 
was  rolled  back  to  305  days,  and  the  production  was  adjusted 
accordingly  for  both  milk  and  fat  to  give  the  total  305-Day 
production.  The  complete  record  to  3I8  days  was  used  in 
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calculating  the  regression,  coefficients,  sums  of  squares,  etc., 
in  the  regression  program. 

Figure  27  gives  the  output  to  the  regression  program. 
Only  the  regression  coefficients  for  milk  and  fat  were  used  in 
the  main  analysis,  however,  the  other  values  were  included  to 
retain  information  which  could  he  used  in  later  studies.  The 
column  headings  should  be  self-explanatory,  as  the  notation  has 
been  used  elsewhere  in  the  study.  The  Y-intercepts  for  milk 
and  fat  are  merely  log,0A  for  milk  and  fat,  and  the  regression 
coefficients  are  the  values  of  C  in  equation  [6]. 


APPENDIX  II 
Over  5O  programs  were  used  to  load  the  data  on  tape , 
list  them,  edit  the  records,  convert  them  into  workable  form, 
etc.  Most  of  these  programs  were  of  such  a  specialized  nature 
that  there  would  be  little  purpose  in  including  them  in  this 
study.   They  were  written  for  the  IBM  1401/1311  system  at  the 
DHI  Computing  Service,  the  IBM  7OAO  at  Brigham  Young  University, 
and  the  IBM  70W7O94  at  the  University  of  Washington.   Several 
of  the  programs  used  are  in  general  use,  however.  The  program 
used  to  test  for  "goodness  of  fit"  of  the  transformation  to 
equation  [6]  was  STAT  03  from  the  Brigham  Young  University 
Computer  Research  Center.   The  generalized  sorting  routine 
from  this  center  was  used  on  all  tape  sorts.   The  program 
BMD02V  of  the  BMD  series  [Dixon  (1965)]  was  used  for  all  the 
factorial  analyses  included  in  the  study.   This  program  was 
run  on  the  IBM  7094  at  the  University  of  Washington.   BMD09S 
was  used  to  convert  to  logarithms  and  antilogs  to  transform 
the  regression  coefficients  into  the  form  in  which  the  factors 
are  presented  by  most  workers  [see  the  introduction]  for  com- 
parisons.  The  regression  program  giving  Figure  2?  as  output 
is  also  quite  specialized,  but  very  simple,  and  therefore  was 
also  not  included. 
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ABSTRACT 

A  study  was  made  on  records  of  18,541  cows  from  North 
Dakota,  Iowa,  Nebraska,  Montana,  Idaho,  Wyoming,  New  Mexico, 
Arizona,  Colorado,  Utah,  and  Nevada,  records  processed  on  a    " 
monthly  basis  by  the  DHI  Computing  Service,  Provo,  Utah. 
Eight  variables  were  studied  to  determine  their  effect  on  the 
shape  of  the  lactation  estimation  curve,  i.e.,  the  curve 
representing  the  factors  used  to  estimate  total  305-day  from 
cumulative  production  to  date.   The  variables,  i.e.,  geographic 
area,  per  cent  fat,  lactation  number,  season,  previous  days 
dry,  days  open,  level  of  milk  production,  and  level  of  fat 
production,  were  broken  into  five  analyses  and  studied  in 
terms  of  their  effect  on  both  the  milk  and  fat  lactation 
estimation  curves. 

The  ratios  of  total  to  part  milk  and  fat  production  were 
computed  for  every  stage  of  the  lactation  in,  terms  of  days  in 
milk,  where  days  in  milk  varied  from  one  to  3O5  days,  for  each 
cow  in  the  study.   The  date  for  each  cow  were  transformed, 
using  a  log-log  transformation,  to  give  one  regression  coef- 
ficient characterizing  the  entire  lactation  estimation  curve. 
In  the  five  factorial  analyses  for  milk  and  fat,  the  common 
regression  coefficient  for  each  cell  was  considered  the 
observation  to  be  analyzed,  with  the  exception  of  the  analyses 
for  level  of  production,  where  the  individual  records  were 
considered  the  observations. 


The  analysis  of  variance,  the  analysis  of  components  of 
variance,  the  analysis  of  means,  and  the  analysis  of  inter- 
actions were  used  in  the  study. 

The  analysis  of  variance  for  the  five  analyses  indicated 
that  Per  Cent  Fat,  lactation  Number,  Season,  and  Level  of  Milk 
and  Fat  Production  were  significant  for  hoth  milk  and  fat 
coefficients.  Days  Open  and  Days  Previous  Dry  were  found  to 
he  not  significant.  Per  Cent  Fat  was  found  to  interact  with 
Season  in  almost  every  case,  in  one  case,  Lactation  Number 
interacted  with  Geographic  Area,  and  in  two  instances,  Per'  Cant 
Fat  interacted  with  Geographic  Area. 

The  analysis  of  the  components  of  variance  indicated  that 
in  all  cases,  Lactation  Number  contributed  the  largest  percent- 
age of  the  variance.  Per  Cent  Fat  and  Season  contributed  more 
to  the  total  variance  for  milk  coefficients  than  for  fat 
coefficients. 

The  multiple  comparisons  of  means  using  Tukey's  hsd  test 
indicated  a  different  structure  for  the  milk  coefficients  than 
for  the  fat  coefficients,  although  these  results,  in  terms  of 
which  levels  could  be  combined  for  each  of  the  variables,  could 
not  always  be  taken  as  representative  because  of  the  presence 
of  significant  interaction. 

The  analysis  of  the  interaction  effects  found  to  be  sig- 
nificant indicated  that  Per  Cent  Fat  by  Season  interactions 
were  quite  consistent  and  important  relationships  both  with 
respect  to  milk  and  to  fat.   It  was  suggested  that  both  Breed 


and  level  of  Average  Per  Cent  Fat  played  roles  in  determining 
both  the  interaction  and  the  per  cent  fat  effect,  because  of 
the  definite  jumps  in  the  graphs  of  the  data  for  the  high  fat 
level  cows. 

Because  of  the  effects  of  the  interaction  in  most  of  the 
relationships,  the  factors  were  given  in  terms  of  the  cell  common 
regression  coefficients  for  each  of  the  original  cells  of  the 
analysis,  with  the  idea  that  further  research  needed  to  he  done 
to  separate  and  study  the  Breed  by  Per  Cent  Fat  by  Season 
effects. 


